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Abstract. We prove that the two-dimensional Schrodinger operator with a potential having the 
symmetry of a honeycomb structure has dispersion surfaces with conical singularities {Dirac points) 
at the vertices of its Brillouin zone. No assumptions are made on the size of the potential. We also 
prove the robustness of such conical singularities to a restricted class of perturbing potentials, which 

■ break the honeycomb lattice symmetry. General small perturbations of potentials with Dirac points 
do not have Dirac points; their dispersion surfaces are smooth. The presence of Dirac points in 

' honeycomb structures is associated with many novel electronic and optical properties of materials. 

(N" 

X) 

1. Introduction and Outline. In this article we study tiie spectral properties 
of the Schrodinger operator Hy = —A + V(x), x e M^, where the potential, V, is 
periodic and has honeycomb structure symmetry. For general periodic potentials the 
spectrum of Hy, considered as an operator on L^(R^), is the union of closed intervals 
of continuous spectrum called the spectral bands. Associated with each spectral band 
are a band dispersion function, /i(k), and Floquet-Bloch states, u(x; k) = p(x; k)e'''", 

■ where i/u(x;k) = /i(k)M(x;k) and p(x; k) is periodic with the periodicity of V{x). 
, The quasi- momentum, k, varies over B, the first Brillouin zone [8]. Therefore, the 

time-dependent Schrodinger equation has solutions of the form e*('^ ''~'^('^)*) p(x;k). 
Furthermore, any finite energy solution of the initial value problem for the time- 
dependent Schrodinger equation is a continuum weighted superposition, an integral 
dk, over such states. Thus, the time-dynamics are strongly influenced by the character 
of ^(k) on the spectral support of the initial data. 
, We investigate the properties of /i(k) in the case where V = Vh is a, honeycomb 

OO ■ lattice potential, i.e. Vh is periodic with respect to a particular lattice, Ah, and has 

CO , honeycomb structure symmetry; see Definition 12.11 There has been intense interest 

fSJ ■ within the fundamental and applied physics communities in such structures; see, for 

, example, the survey articles [12l[T4]. Graphene, a single atomic layer of carbon atoms, 

is a two-dimensional structure with carbon atoms located at the sites of a honeycomb 
structure. Most remarkable is that the associated dispersion surfaces are observed to 
have conical singularities at the vertices of B^, which in this case is a regular hexagon. 
That is, locally about any such quasi- momentum vertex, k K,, one has 



C^: /i(k)-MK.) * + |A„|.|k-K.| , (1.1) 

for some complex constant Aj 0. A consequence is that for wave-packet initial 
conditions with spectral components which are concentrated near these vertices, the 
effective evolution equation governing the wave-packet envelope is the two-dimensional 
Dirac wave equation, the equation of evolution for massless relativistic fermions [T2J[T]. 
Hence, these special vertex quasi-momenta associated with the hexagonal lattice are 
often called Dirac points. In contrast, wave-packets concentrated at spectral band 
edges, bordering a spectral gap where the dispersion relation is typically quadratic, be- 
have as massive non-relativistic particles; the effective wave-packet envelope equation 
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is the Schrodinger equation with inverse effective mass related to the local curvature of 
the band dispersion relation at the band edge. The presence of Dirac points has many 
physical implications with great potential for technological applications [50]. Refrac- 
tive index profiles with honeycomb lattice symmetry and their applications are also 
considered in the context of electro-magnetics [5l [19] . Also, linear and nonlinear prop- 
agation of light in a two-dimensional refractive index profile with honeycomb lattice 
symmetry, generated via the interference pattern of plane waves incident on a pho- 
torefractive crystal, has been investigated in [HI [3] . In such structures, wave-packets 
of light with spectral components concentrated near Dirac points, evolve diffractively 
(rather than dispersively) with increasing propagation distance into the crystal. 

Previous mathematical analyses of such honeycomb lattice structures are based 
upon extreme limit models: 

1. the tight-binding / infinite contrast limit (see, for example, [TB] [T^l E] ) in 
which the potential is taken to be concentrated at lattice points or edges 
of a graph; in this limit, the dispersion relation has an explicit analytical 
expression, or 

2. the weak-potential limit, where formal perturbation theory is used [5l[l]. 

The goal of the present paper is to provide a rigorous construction of conical 
singularities (Dirac points) for essentially any potential with a honeycomb structure. 
No assumptions on smallness or largeness of the potential are made. More precisely, 
consider the Schrodinger operator 

H'^^^ = -A + eVh (ereal) (1.2) 

where 14 (x) denotes a honeycomb lattice potential. These potentials are real- valued, 
smooth, A/i- periodic and, with respect to some origin of coordinates, inversion sym- 
metric (x — x) and invariant under a 27r/3- rotation {TZ- invariance); see Def. 12.11 
We also make a simple, explicit genericity assumption on Vh{x); see equation (|5.2p . 

Our main results are: 

1. Theorem 15. Ji which states that for fixed honeycomb lattice potential Vh, the 
dispersion surface of H^^^ has conical singularities at each vertex of the hexag- 
onal Brillouin zone, except possibly for e in a countable and closed set, C. We 
do not know whether exceptional non-zero e can occur, i.e. whether the above 
countable closed set can be taken to be {0}. However our proof excludes ex- 
ceptional e from (— eo, £o)\{0}, for some Eq > 0. Moreover, for small e these 
conical singularities occur either as intersections between the first and second 
band dispersion surfaces or between the second and third dispersion surfaces. 
As £ increases, there continue to be such conical intersections of dispersion 
surfaces, but we do not control which dispersion surfaces intersect. 

2. Theorem \9.1\ which states that the conical singularities of the dispersion sur- 
face of H^'^^ fore ^ C, are robust in the following sense: Let W be Ah- periodic 
and inversion-symmetric, but not necessarily TZ- invariant. Then, for all suf- 
ficiently small real rj, the operator H{r]) = H^'^'^ + riW has a dispersion surface 
with a conical-type singularity. Furthermore, these conical singularities will 
typically not occur at the vertices of the Brillouin zone, Bt; see also the nu- 
merical results in In Remark 19.21 we show instability of Dirac points to 
certain perturbations, e.g. perturbations W which are Kh- periodic but not 
inversion-symmetric. The dispersion surface is smooth in this case. 
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The paper is structured as follows. In section [5] we briefly outline the spectral 
theory of general periodic potentials. We then introduce A^, the particular lattice 
(Bravais lattice) used to generate a honeycomb structure or "honeycomb lattice" , 
the union of two interpenetrating triangular lattices. Section [2] concludes with im- 
plications for Fourier analysis in this setting. Section [3] contains a discussion of the 
spectrum of the Laplacian on L^, the subspace of satisfying pseudo-periodic bound- 
ary conditions with quasi- momentum k e Sh, the Brillouin zone. We observe that 
degenerate eigenvalues of multiplicity three occur at the vertices of Bh- In section 
m we state and prove Theorem 14.11 which reduces the construction of conical singu- 
larities of the dispersion surface at the vertices of Bh to establishing the existence 
of two-dimensional TZ— invariant eigenspaces of for quasi-momenta at the ver- 
tices of Bh- In section [5] we give a precise statement of our main result, Theorem 
15.11 on conical singularities of dispersion surfaces at the vertices of Bh- In section [S] 
we prove for all e sufficiently small and non-zero, by a Lyapunov-Schmidt reduction, 
that the degenerate, multiplicity three eigenvalue of the Laplacian splits into a multi- 
plicity two eigenvalue and a multiplicity one eigenvalue, with associated TZ- invariant 
eigenspaces. In order to continue this result to e large we introduce, in section [Tj 
a globally-defined analytic function, £(^,e), whose zeros, counting multiplicity, are 
the eigenvalues of H'-'^\ Eigenvalues occur where an operator / + C{iJ,,e), C{ii,e) 
compact, is singular. Since C(/i,£) is not trace-class but is Hilbert-Schmidt, we work 
with e) = det2(/ + C(/i, e)), a renormalized determinant. In section[8l f (/i, e) and 
A| (see (11.11) ) are studied using techniques of complex function theory to establish the 
existence of Dirac points for arbitrary real values of e, except possibly for a countable 
closed subset of M. In section IH] we prove Theorem 19.11 which gives conditions for 
the local persistence of the conical singularities. Remark 19.21 discusses perturbations 
which break the conical singularity and for which the dispersion surface is smooth. 
Appendix |X] contains a counterexample, illustrating the topological obstruction dis- 
cussed in section [831 
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like to thank Z.H. Musslimani for stimulating discussions early in this work. We are 
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1.1. Notation. 

1. z e C =^ z denotes the complex conjugate of Z- 

2. A, a. d X d matrix =^ A* is its transpose and A* is its conjugate-transpose. 

3. <j> = VTTW- 

4. K"" = K™i'™2 = K + mk = K + miki + TO2k2. 
K,ki and k2 are defined in section [221 

5. Vk = e-^kocy^gzk.x = Vx + ik, Ak = Vk • Vk. 

6. X, y e C", <x, y> = x • y, x • y = xiyi + ■- - + XnVn- 

7. For q = (qi_, 92) e Z^, qk = q^ki + q2k2 - 

8. </,3> = //ff 

9. = = : Sj,^. <j>2iejP < ^} 

2. Periodic Potentials and Honeycomb Lattice Potentials. We begin this 
section with a review of Floquet-Bloch theory of periodic potentials [1] , [TU] , [TB] . We 
then turn to the definition of honeycomb structures and their Fourier analysis. 



4 Honeycomb Lattice Potentials and Dirac Points 

2.1. Floquet-Bloch Theory. Let {vi, V2} be a linearly independent set in M?. 
Consider the lattice 

A = {mivi + m2V2 : mi, m2 e Z } = Zvi © Zv2, (2.1) 

The fundamental period cell is denoted 

n = { 6*1 vi + 6I2V2 : ^ 6lj ^ 1, j = 1, 2 } . (2.2) 

Denote by L^^^ ^ = L^(R^/A), the space of Lf^^ functions which are periodic with 
the respect to the lattice A, or cquivalcntly functions in on the torus M^/A = T^: 

/ G Lper.A if ^^d Only if /(x + v) = /(x), for x 6 R^, v e A . 

More generally, we consider functions satisfying a pseudo-periodic boundary condition: 

/ e Ll f^ if and only if /(x + v) = f{x)e^-^, for x e R^ v e A. (2.3) 

We shall suppress the dependence on the period-lattice. A, and write L^, if the choice 
of lattice is clear from context. For / and g in L^j^, fg is locally integrable and A- 
periodic and we define their inner product by: 

(f,g) = / 7R5(x) dx. (2.4) 
Jn 

In a standard way, one can introduce the Sobolev spaces ^. 
The dual lattice. A*, is defined to be 

A* = {miki + m2k2 : mi, m2 e Z} = Zki © Zk2 , (2.5) 

where ki and k2 are dual lattice vectors, satisfying the relations: 

ki • Vj = 2Tr6ij . 

If / e ipgr then / can be expanded in a Fourier series with Fourier coefficients 

/ = {/m}meZ2: 

/(X) = 2 = E /-i.™^ ei(mxki+,n.k.).x ^ (2.6) 

^ ^ ^ e-™''-^ /(y) dy = |1 ^ g-.(™,ki+,n.k.).y _^(y) (2.7) 

Let V^(x) denote a real- valued potential which is periodic with respect to A, i.e. 

y(x-F v) = F(x), for xe R^ veA 

Throughout this paper we shall also assume the potential, V{x) , under consideration 
is C". Thus, 

VeC^iR^/A). (2.8) 

We expect that this smoothness assumption can be relaxed considerably without much 
extra work. 
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For each k e we consider the Floquet-Bloch eigenvalue problem 

Hv 0(x,k) = ^(k) 0(x,k), xe (2.9) 

0(x + v,k) = e*''-^ 0(x,k), veA, (2.10) 

where 

Hv = -A + y(x) . (2.11) 

An L^- solution of (|2.9I) - (I2.10I) is called a Floquet-Bloch state. 

Since the eigenvalue problem (|2.9p - (|2.10p is invariant under the change k i-^ k + k, 
where k e A*, the dual period lattice, the eigenvalues and eigenfunctions of (|2.9p - 
(|2.10p can be regarded as A*— periodic functions of k, or functions on = M^/A*. 
Therefore, it suffices to restrict our attention to k varying over any primitive cell. It 
is standard to work with the first Brillouin zone, B, the closure of the set of points 
k e R^, which are closer to the origin than to any other lattice point. 

An alternative formulation is obtained as follows. For every k e ;B we set 

0(x; k) = e*'' ''p(x; k) (2.12) 

Then p(x; k) satisfies the periodic elliptic boundary value problem: 

ffy(k)p(x,k) = ^(k)p(x,k), xeM^ (2.13) 

p(x + v;k) =p(x,k), veA, (2.14) 

where 

Fy(k) = -(V + ik)Vl^(x) = -Ak + F(x) . (2.15) 

The eigenvalue problem (|2.9p - (l2.10l) . or equivalently (I2.13p - (|2.14p . has a discrete spec- 
trum: 

Aii(k) < Ai2(k) < Ai3(k) < ... (2.16) 

with cigenpairs pf,(x; k), /ib(k) : = 1, 2, 3, . . . . The set {p6(x; k)}fc^i can be taken 
to be a complete orthonormal set in Lp(,^(R^/A). 

The functions ^fc(k) are called band dispersion functions. Some general results 
on their regularity appear in [2]. As k varies over B, /Ltb(k) sweeps out a closed real 
interval. The spectrum of —A + ^(x) in L^(R^) is the union of these closed intervals: 

spec(i/y) = y spec(i/y(k)) . (2.17) 

Moreover, the set lJ6>i Ukes{'^'>(-''-' '^''(^j ^) = e'''"Pb(x; k), suitably normalized, 
is complete in L^(R^): 

/e l2(m2) ^ ^(^) = E / <4(-,k),/>0b(x,k) dk , 
where the sum converges in the norm. 



6 Honeycomb Lattice Potentials and Dirac Points 

2.2. The period lattice, A/j , and its dual, A,*. 

Consider A/i = Zvi ©Zv2, the lattice generated by the basis vectors: 

vi = a I 1 V2 = a I , a > 0. (2.18) 

The dual lattice A^ = Zki © Zk2 is spanned by the dual basis vectors: 



1 \ /I 

2 \ / 2 



47r 



, , k2 = g , 9 - ^ , (2.19) 

V3 / \ _V3 / aV3 



2 / \ 2 



where 



• ve, = 27r(5«- , (2.20) 

|vi| = |v2|=a, Vi-V2 = y , (2.21) 

|ki| = |k2| = (Z, ki.k2 = -i<z2 . (2.22) 

The Brillouin zone, Bh, is a hexagon in M^; see figure \2A\ Denote by K and K' 
the vertices of Bh given by: 

K = i (ki - k2) , K' = -K = i (k2 - ki) . (2.23) 

All six vertices of Bh can be generated by application of the rotation matrix, R, which 
rotates a vector in clockwise by 27r/3. R is given by 



1 


V3 


2 


2 


V3 


1 


2 


2 



i? = (2.24) 



and the vertices of Bh fall into to groups, generated by action of i? on K and K': 
K type - points : K, RK = K + k2, R^K = K - ki 
K' type - points : K', RK' = K' - k2, R^K' = K' + ki . (2.25) 

Remark 2.1 (Symmetry Reduction). Let ( 0(x, ; k), /i(k) ) denote a Floquet- 
Bloch eigenpair for the eigenvalue problem (|2.9p - (|2.10p with quasi-momentum k. Since 
V is real, ( 0(x; k) = (/)(x; k), /i(k) ) is a Floquet-Bloch eigenpair for the eigenvalue 
problem with quasi-momentum — k. Recall the relations (|2.25p and the A^- periodicity 
of: k I— > /i(k), are k i— > 0(x;k). /t follows that the local character of the dispersion 
surfaces in a neighborhood of any vertex of Bh is determined by its character about 
any other vertex of Bh ■ 

In our computations using Fourier series, we shall frequently make use of the 
following relations: 

i?ki=k2, i? k2 = - (ki + k2) , i? (ki + k2) = -ki (2.26) 

RQ= T RC = f C, where (2.27) 

I f I \ 2^ 1 . .V3 3 



C=^^, ), r^e^^-- + ^-, ^ 1. (2.28) 
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Moreover, R* maps the period lattice A/j to itself and, in particular, 

i?*Vi = -V2, i?*V2 = Vi - V2 (2.29) 

2.3. Honeycomb lattice potentials. 

For any function /, defined on M^, introduce 

7e[/](x) = /(i?*x), (2.30) 

where R is the 2x2 rotation matrix displayed in (|2.24[) . 
Definition 2.1. [Honeycomb lattice potentials] 

Let V be real-valued and V e C"^(R^). V is a honeycomb lattice potential if there 
exists Xo e such that V(x.) = ^(x — Xq) has the following properties: 

1. V is A/i— periodic, i.e. V{x + v) = V^(x) for all x e and v e A^. 

2. V is even or inversion-symmetric, i.e. V{—x) = T^(x). 

3. V is TZ- invariant, i.e. 

TZ[V]{x) = V^(i?*x) = f(x), 

where, R* is the counter-clockwise rotation matrix by 27r/3, i.e. R* = R^^, 
where R is given by (j2.24p . 

Thus, a honeycomb lattice potential is smooth, A^- periodic and, with respect to 
some origin of coordinates, both inversion .symmetric and TZ- invariant. 

Remark 2.2. As the spectral properties are independent of translation of the 
potential we shall assume in the proofs, without any loss of generality, that xq = 0. 

Remark 2.3. A consequence of a honeycomb lattice potential being real-valued 
and even is that if ((/)(x; k), /i) is an eigenpair with quasimomentum k of the Floquet- 



Bloch eigenvalue problem, then |^0(— x; k), /ij is also an eigenpair with quasimomen- 
tum k. 

Remark 2.4. We present two constructions of honeycomb lattice potentials. 
Example 1: "Atomic" honeycomb lattice potentials: Start with the two points 

A =(0,0), and B = a^--^,0^ , (2.31) 

which lie within the unit period cell of A^; see (j2.18[) . Define the triangular 
lattices of A- type and B- type points: 

Aa = A + A/j and Ab = B + A,^ (2.32) 

We define the honeycomb structure, H, to be the union of these two triangular 
lattices: 

H = Aa u Ab ; (2.33) 

see Figure \2.1\ Let Vq be a smooth, radial and rapidly decreasing function, 
which we think of as an "atomic potential" . Then, 



V{^)= K)(x-a), 
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is a potential associated with "atoms" at each site of the honeycomb structure 
H. Moreover, V^(x) is a honeycomb lattice potential in the sense of Definition 
\K^with xo = B. 

Example 2: Optical honeycomb lattice potentials: The electric field envelope of a 
nearly monochromatic beam of light propagating through a dielectric medium 
with two-dimensional refractive index profile satisfies a linear Schrodinger 
equation idzip = — Aj. yi/i + V{x, y)il) = 0. Here, z denotes the direction 
of propagation of the beam and {x, y) the transverse directions. Honeycomb 
lattice potentials have been generated taking advantage of nonlinear optical 
phenomena. It was demonstrated in 115] that a honeycomb lattice potential 
(a honeycomb "photonic lattice"), V{x,y), can generated through an optical 
induction technique based on the interference of three plane wave beams of 
light within a photorefractive crystal, exhibiting the defocusing (nonlinear) 
optical Kerr effect. The refractive-index variations are governed by a potential 
of the approximate form: 



y(x) * Vo ( cos(ki • x) + cos(k2 • x) + cos ((ki + k2) • x)) 



(2.34) 



It is straightforward to check, in view of (|2.26p . that a potential of this type 
is a honeycomb lattice potential in the sense of Definition \2.1\ with xq = 0. 
In fact, in Proposition \2.3\ below we assert that with respect to some origin 
of coordinates, any honeycomb lattice potential can be expressed as a Fourier 
series of terms of this type. 





0.5 1 



Fig. 2.1. Left plot: Part of the honeycomb structure, H. H is the union of two sub-lattices 
Aa = a + Afi (blue) and Ab = B + A^ (green). The lattice vectors {vi,V2} generate A/j. See 
Remark \2.4\ Right plot: Brillouin zone, B^, and dual basis {ki,k2}. K and K' are labeled. Other 
vertices of Bh obtained via application of R, rotation by 27r/3; see equation 112.251 1 . 



The following proposition plays a key role. It states that at distinguished points 
in k— space, namely the K and K' type points, Hv with quasi-momentum depen- 
dent boundary conditions (|2.10|) or equivalently, iJy(k), with periodic boundary 
conditions, has an extra rotational invariance property. 

Proposition 2.2. Assume V is a honeycomb lattice potential, as in Definition 
\2.1l Assume K* is a point o/K or K' type; see (I2.25|) . Then, H and TZ map a dense 
subspace of to itself. Furthermore, restricted to this dense subspace of , the 
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commutator [H, TZ\ = H TZ — TZ H vanishes. In particular, if (/'(x, k) is a solution of 
the Floquet-Bloch eigenvalue problem (I2.9I) - (|2.10I) luith k = K*, then ??.[(/)(•, k)](x) is 
also a solution of ((?!^ - (|2.10p with k = K*. 

Proof. Take as a dense subspace , the space of functions satisfying /(x + 
v) = e*^*-^/(x) for all x e and v s'a^. Clearly, H maps C£ to itself. Define 
0fl;(x) = ??.[(/>(•, K,)](x) = (?!)(i?*x, K*). Without loss of generality, assume K» = K. 
By (12211), if V e A,j then iT!*v e K,,. We have 

^^(x + v) = 0(i?*x+i?*v,K) = e'^-^*"(/)(i?*x,K) 

= e*-^^-^ 0(i?*x,K) = e*(^+''^>^ (/)(i?*x,K) 
= e*^-^ 0(i?*x,K) = e'^-^ 0h(x). 

Thus, we have 7?. maps to itself. 

Next note that by invariance of the Laplacian under rotations, — Ax</!)i?(x) = 
— Ay(/)(y, K«)|y^^,|,^ . Furthermore, by TZ— invariance of l^(x), that V{'x)(f)ii{yL) = 

y(i?*x)0(i?*x,K*) = V"(y)<?!)(y,K*)|y^^*x • Therefore, [i^,7^] vanishes on on Cg^. 
In particular, we have that 

i7y0(x, K«) = /x(/)(x, K*) =^ Hv(I>r{-^) = f^(f>R{:>i-) ■ 

This completes the proof of the proposition. □ 

We conclude this section with a discussion of the Fourier representation of hon- 
eycomb lattice potentials in the sense of Definition 12. II Let ^(x) be such a potential 
with Fourier series: 

meZ2 (mi,m2)eZ2 

Since V(x) = 7^[F](x), we have 

T/Zp* \ _ T/ i(mii?,ki+m2-Rk2)-x _ t/ — m2)ki + (mi — m2)k2)-x 

(mi,m2) (mi,m2) 

Therefore, V;„i,m2 = T^-m2.rni-m2- Similarly, V^(x) = 7?.^[F](x) implies that F„ii,m2 = 
Vm2-mi,-mi- lutroducc the mapping i? : — > acting on the indices of the Fourier 
coefficients of V: 

i?(mi,TO2) = (— TO2,'Tii — m2) and therefore 

_R^(mi, TO2) = (m2 — mi,— mi), and R^{mi,m2) = (mi,m2) . (2.35) 
Then we have 

= Vf,^ = y^2^ (2.36) 

Note that ^0 = and that is the unique element of the kernel of R. Furthermore, 
any m 7^ lies on an R— orbit of length exactly three. Indeed, 

m = i?m <-> (mi, m2) = (— m2, mi — m2) =^ mi = m2 = and 
m = i?^m <-> (7711, m2) = (— mi + m2, — mi) =^ mi = m2 = . 
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Suppose m and n are non-zero. We say that m ~ n if m and n lie on the same 3— 
cycle. The relation ~ is an equivalence relation, which partitions Z\{0} into equiva- 
lence classes, (Z\{0}) / ~. Let S denote a set consisting of exactly one representative 
from each equivalence class. We now have the following characterization of Fourier 
series of honeycomb lattice potentials: 

Proposition 2.3. Let V(x) denote a honeycomb lattice potential. Then, 
V{x) = V{0)+ l^cos(mk-x) + cos((:^m)k • x) + cos((7e2m)k • x) , 

me 5 

(2.37) 

where Vm are real and R is defined in (j2.35p . 
Proof of Proposition \2.3[ From (|2.36p we have 

y(x) = V{0) + XI ^™ l^gimk-x _^ gi(J^m)k.x _^ g.(fl'm)k.x^ ^2.38) 

The relation V{x.) = (F(x) + V^(-x))/2 and ([Qg]) imply 
V{x) = V{0) + X (cos(mk • x) + cos ({Rm)k ■ icj + cos (^{R^m)k ■ x^ j 

(2.39) 

Moreover, since V is real and even, Vm is real for me 1?. This completes the proof. 
2.4. Fourier analysis in . 

We characterize the Fourier series of functions (p 6 -^k > ^- ^- functions 0(x; K« ) , 
satisfying the quasiperiodic boundary condition: 

0(x + v) = e*^*-" 0(x) (2.40) 

The discussion is analogous to that preceding Proposition 12.31 

If p.40p holds then 0(x) = e*^*'^ p(^)i where p(x) is A— periodic. It follows 
that (/) has a Fourier representation: 

0(x) = e^^*-" 2 c(mi,m2) e^(™i*'i+™^*'^>", (2.41) 

(mi ,7712)6^2 

which we rewrite as 

^(x) = Yi c(mi,m2) e*(^*+™i*'i+™^''^^-'' 

2 c(TOi,m2) e*^""'"'-" = ;X c(m)e''K"-^, (2.42) 

where K™ = + miki + m2k2. 
Usually, we denote by c^ijn.) or c(m; 0) t/ie Fourier coefficients, as in p. 421) . 0/ 

Note that the transformation TZ, defined in p.30p . is unitary on and so its 
eigenvalues lie on the unit circle in C. Furthermore, if TZcj) = Xcj) and (f> ¥= 0, then 
since TZ^ = Id, (f) = TZ^cj) = X^(f), we have A'^ = 1. Therefore A e {l,r, f}, where 
r = exp(27ri/3). 
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We are interested in the general Fourier expansion of functions in each of the 
eigenspaces of TZ: 

= {/ e ^K. ■T^f = f} (2.43) 
iK.,r = {/eL^. :7e/ = r/}, (2.44) 

^K... = {/ e • ^/ = ^/} (2-45) 

Since TZ is unitary these subspaces are pairwise orthogonal. 

Fix, without loss of generality, K« = K. We first consider the action of TZ on 
general </> e L^^. Applying TZ to (f>, given by (|2.42p . we obtain: 

TZm^)= ^ c4mi,m2) e^^'"-^*^ 

(mi ,7712)^^2 
(mi ,m2)GZ^ 

2 C0(mi,m2) e*(K+(-"^)''i+('"i-™^+i)'^^)-", 

(mi,m2)eZ2 

since 

RK"" = i^K™!^"^ = K + (-m2)ki + (mi - ma + l)k2 = K-m2,n.i-m^+i ^ (2.46) 
Thus, 

Cizcj,{-~rn2,'rni — m2 + 1) = C0(mi,m2), or equivalently 

C7?,0(mi, TO2) = C0(m2 - mi - 1, -mi) . (2.47) 

Similarly, by a second application of TZ, and using the relation 

J^2j^mi,m2 _ J^m2— mi — 1, — mi (2 48^ 

we have 

C7j20(m2 — mi — 1, —mi) = c^{mi, 1112), or equivalently 

C7?,20(mi, ma) = c^{-m2, mi - m2 + 1) (2.49) 

Finally, since = I, c-jis ^{1711,1712) = €^{1711, 1712). 

TZ acting in induces a decomposition of 7? into orbits of length three: 

(mi,m2)^i— > (— m2, mi — m2 + 1) ^ I— > (ma — mi — 1, — mi) ^ i-^ (mi,m2) (2.50) 

For convenience we shall abuse notation and write 

TZm = 7?,(mi, ma) = (—ma, mi — ma + 1) 
TZ^m = 7?,^(mi, ma) = (ma — mi — 1, —mi), 
TZ^m = Id (mi, ma) = (mi, ma) (2-51) 

Using the notation (|2.5ip . relations (|2.47l) . (|2.49l) and (|2.5ip can be expressed as: 

CTC0(m) = C4,(TZ^ui) = C0(ma - mi - 1, -mi) 

C7j20(m) = c^{TZm) = C0(-ma, mi - ma + 1) (2.52) 
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Furthermore, by ((^^ and ((^^ 

i?K'" = K^'" and R^K"" = K'^'"'. (2.53) 

Each point in Hes on an orbit of TZ of precisely length 3, a 3-cycle. To see 
this, note that by (|2.51l) TZ^m = m for all m e Z^. So we need only check that 
there are no solutions to either TZm = m or to TZ^m = m. First, suppose TZm = m. 
Then, TZ^m = m as well. So, on the one hand the centroid of m, TZm and TZ^m is 
equal to m e Z^. On the other hand, by (|2.5ip their centroid is (—1/3, 1/3) ^ Z^, a 
contradiction. Therefore, there are no Z^ solutions of TZm = m. Now if TZ^m = m, 
then applying TZ to this relation yields m = T^.'^m = TZm, and we're back to the 
previous case. 

We shall say that two points in 7?, m and n are equivalent, m ~ n, if they lie 
on the same 3-cycle of TZ. We identify all equivalent points by introducing the set of 
equivalence classes, 1? j ~ . 

Definition 2.4. We denote by S a set consisting of exactly one representative of 
each equivalence class in 1? j ~ . For example, {(0, 0), (0, 1), (— 1, 0)} e 1? / ~, from 
which we choose (0, 1) as its representative in S. 

Using the relations (|2.5ip , we can express the Fourier series of an arbitrary e 
as a sum over 3-cycles of TZ: 

0(x) = 2 c^(mi,m2) e'(K+"i'^i+"'^''^>'' 

(mi ,m2)eZ^ 

= (c0(m) e*^"-^ + c47em)e*«^""-^ + C4,{TZ'm)e'''"'^"'-^ ) , (2.54) 

me 5 

where TZ^m, j = 1,2 is given in (I2.5ip . 

We now turn to the Fourier representation of elements of the subspaces ^ , ^ 
and i^^. 

Proposition 2.5. Let 4> e Lk- 

TZ(j) = (j) ^ C0(m) = c^{TZm) = Cfj^iTZ^m) 
TZ(j) = T (f) c^{TZ^m) = Tc^{m) and c^{TZm) 
TZ(f) = f (f) c^{TZ^m) = fc^{m) and c^{TZm) 

Moreover, 

TZ^(j) = f (p <=> CK20(m) = C0(7^m) = fc^{m), (2.58) 

where TZm and TZ^m are defined in (|2.5ip . 

Proof Assume TZcp = acj). Then, CTi<j,{m) = ac^{m). By (|2.52p CTi^{m) = 
c^{TZ'^m) and therefore 

C0(7^2m) = cr C0(m) (2.59) 

Furthermore, 7?.^0 = cr'^cj) and therefore CTi2^{m) = (T^C0(m). By p.52p CTz2^[m) = 
c,f,{TZm) and therefore 

C0(7^m) = cr^ C0(m) (2.60) 

We next apply relations (I2.59P and (I2.60p to the cases: a = 1,t, f. Let a = t. 
Then, TZ(f> = rcj) implies C4,{TZ'^m) = r C0(m), by ((^3^)) . Also, by ((^^ we have 



(2.55) 

= fc^(m) (2.56) 
= TcJm). (2.57) 
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c^{TZm.) = C0(m) = r c^{ni). This proves (|2.56p . The cases a = 1,t are similar. 
□ 

Proposition 12 . 51 can now be used to find a representation of the eigenspaces of TZ. 
We state the result for an arbitrary point, K*, of K or K' type. 



Proposition 2.6. 

1. (p e ^ <^ there exists {c(m)}mes £ l^i^) such that 

^(x) = 2 c(m) (e'^" " + f e'^^" " + re'-"'^"" ) . (2.61) 

meS 

2. 4> e ^ o there exists {c(m)}„ies e P{S) such that 

0(x) = 2 c(m) (e'K""' + re'^K"-^ + fe^^'K^x ^ (2.62) 

me 5 

5. // e ^ is given by 

</.i(x,K.) = 2 c(m) (e*^" " + f e'^^^ '' + r e^^'^^ ") , (2.63) 

me 5 



then (/)2(x, K*) = (j!)i(— x, K*) e - and 
02(x,K,) = 2 c(^i^) (e*^"-" + re^^-K"-" + f e*^'^"'") . (2.64) 

me5 

4- 4> ^ 1 i/jere exists {c(m)}me5 e ^^(i^) sttc/i i/iat 

0(x) = c(m) (e'^" '' + e'-^K"-^ + e'-^'K^ x ^ _ ^2.65) 

me 5 

We summarize the preceding in a result which facilitates the study of H on 
in terms of the action of TZ on invariant subspaces of H . 

Proposition 2.7. Let K* denote a point of K or K' type, R denote the 27r/3 
clockwise rotation matrix (see (|2.24p ) and 7?,[/](x) = f(R*x). Then TZ, acting on 
has eigenvalues 1, r and f = inducing a corresponding orthogonal sum de- 
composition of into eigenspaces: 

^K. =^K.,l©^K.,r©^K.,^ • (2-66) 

The elements of each summand are represented as in Provosition \2.b\ 

Remark 2.5. Since, by Provosition \2.2\ H and TZ commute on , the spectral 
theory of H in can be reduced to its independent study in each of the eigenspaces 
in the orthogonal sum (|2.66p . 



14 Honeycomb Lattice Potentials and Dirac Points 

3. Spectral properties of H^^^ in - Degeneracy at K and K' points . 

Our starting point for the study of H'-'^^ on is the study of i?^") = —A. Consider 
the eigenvalue problem 

Equivalently 0(o)(x;k) = e^^^'V^Hxik), where p^°\-]k) e L^{R'^/Ah): 

F(o)(k)p(o) = -(V + *k)2pW = ^(o)(k)p(o), (3.1) 
p(")(x + v,k) =p(")(x,k), V6 . (3.2) 

The eigenvalue problem p.ip , p.2p has solutions of the form: 

with associated eigenvalues 

^^^°\,m,0^) = |k + 7niki + m2k2|', keB. (3.3) 



Proposition 3.1. Let k = K* denote any vertex of the hexagon Bh (points of 
K orK' type); see Then, 

1. z^^*^^ = |K*p is an eigenvalue of of multiplicity three with corresponding 
three-dimensional eigenspace 

span { e'^* "" , e'^^*-^ , g»i?^K*.x | ^ ^3 4^ 

2. Restricted to each of the TZ— invariant subspaces of 

/las an eigenvalue /i^"-* = |K*p of multiplicity one with eigenspaces: 



span 




■X ^ gji?K, x ^ gifl^K. xj 


c 




span 




■x^ -g.flK.-x^^g.fl^K.. 


''I 


cz L^^^ and 


span 




- + Te*«K,.x^ -g^i?^K,. 


X 





5. /i^'') is the lowest eigenvalue of iJ^^^ in . 

Proof. Without loss of generality, let K* = K. Since R is orthogonal, |K| = 
|i?K| = \R^K\. Therefore, -A^- = |K|2* for * = e'^'^, e*-^^'' and e*'^'^''. It 
follows that = |K|2 is an eigenvalue of multiplicity at least three. To show that 
the multiplicity is exactly three, we seek all m = (mi, 7712) for which |K™p = |Kp. 
Using K™ = K + TOiki + m2k2, we obtain 

nif + m\ + mi — m2 — mim2 = 

which can be zero only if m = (0,0), (0, 1) or (—1,0). In the first instance, K^"'") = 
K. If m = (0,1), then K^^^) = K + ^2 = KK. Finally if m = (-1,0) then 
K(-i,o) = K - ki = (K + ka) - (ki + ka) = i?K + i?k2 = i?(K + ks) = R^K. This 
proves conclusion 1. Proposition 12.61 above, which characterizes the Fourier series 
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of functions in L^^, a = 1,t, r implies conclusion 2. Conclusion 3. holds because 
ml + m2 — r7iiTO2 + mi — ?7i2 5= I for (mi, 7712) e 7? other than (0, 0), (0, 1) and (—1,0). 
□ 

Recall that for each he Bh, the eigenvalues of H^^^ are ordered (|2.16p : 

(k) ^ (k) ^ /.f (k) ^ (k) < . . . (3.5) 

For k = K we have 

|Kp = ^,'(>\K) = (K) = f,f\K) < 4")(K) < . . . (3.6) 
We shall see in section [6] that for small £, the spectrum perturbs to 
either 

(K) = f,i'^ (K) < fi^^ (K) < A^i^) (K) . . . (3.7) 

or 

fi['^ (K) < (K) = (K) < mI^) (K) *S . . . (3.8) 

In either case, the multiplicity three eigenvalue splits into a multiplicity two eigenvalue 
and a simple eigenvalue. The connection between the double eigenvalue and conical 
singularities of the dispersion surface is explained in the next section; see Theorem 
14.11 We shall see from Theorem 15.11 or rather its proof (in section [6]) that for all 
small e, conical singularities occur at all vertices K* of Bh, and that these occur at 
the intersection point of the first and second band dispersion surfaces in the case of 
p.7p . and at the intersection of the second and third bands in the case of p.8p . As 
s increases, we continue to have such conical intersections of dispersion surfaces, but 
we do not control which dispersion surfaces intersect. 

4. Multiplicity two eigenvalues of H^^^ and conical singularities. Let 

K* a point of K or K' type. In this section we show that if H acting in has a 
dimension two eigenspace E^- ©E^^, where E^- and Et= are dimension one subspaces of 
r ^^'^ T' respectively, then the dispersion surface is conical in a neighborhood 
of K.. Recall that we assume V e C""(M7A/,). 



Theorem 4.1. Let H = —A + V , where V^(x) is a honeycomb lattice potential 
in the sense of Definition \2.1\ Let K, denote any vertex of the Brillouin zone, Bh- 
Assume further that 

(hl.r) H has an ^- eigenvalue, /io = m(K*), of multiplicity one, with corre- 
sponding eigenvector $i(x) = $i(x;K,), normalized to have L^(Q,) norm 
equal to one. 

(hl.f) H has an L^ ^- eigenvalue, /io = m(K«), of multiplicity one, with corre- 
sponding eigenvector $2(x) = <I'i(— x), 
(h2) /io = m(K*) is not an eigenvalue of H on L^^^ 1. 
(hS) the following nondegeneracy condition holds: 

Aj = 3 X area(17) x ^ c(m;$i)2 ( M ■ ^ 0, (4.1) 
where {c(m; $i)}m€5 Fourier coefficients o/$i, as defined in Proposition 
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Then H acting on has a dispersion surface which, in a neighborhood o/ k = K» , 
is conical. That is, for k — K» near 0, there are two distinct branches of eigenvalues 
of the Floquet-Bloch eigenvalue problem with quasi-momentum, k; 

^i+{k) - fi{K.) = +|A„| |k-K.| ( l + i?+(k-K.) ) (4.2) 
M-(k)-MK.) = -|A„| |k-K.| ( l + i?„(k-K.) ), (4.3) 

where E+(k) = 0{\k\) as \k\ — > and E+ are Lipschitz continuous functions in a 
neighborhood ofO. 
Remark 4.1. 

1. Elliptic regularity implies that the eigenfunctions $j, j = 1,2 are in i7^(M^). 
Therefore, X;me5(l + |mp)^|c(m)p < co. We conclude that the sum defining 
Aj converges. 

2. In section\^ we study the case of "weak" or small potentials, i.e. V = eVh 
with £ small. For all e such that < jej < e", where e° is a sufficiently small 
positive number, we will: 

(i) verify the double eigenvalue hypothesis (hi) of Theorem \4.1\ by showing the 
persistence of a double eigenvalue due to intersection of the bands one and two 
in case (|3.7[) or bands two and three in case (|3.8p . (ii) verify hypothesis (h2) 
of Theorem \4-.l\ by showing, via explicit calculation, that the ^ eigenvalue 
of H , differs from the double eigenvalue, and (Hi) verify (h3) by showing: 
|A|p = 16 area(17)2 TT^/a^ + 0{e): see ([QH)) . Theorem\Jl\ then implies 
the existence of a non-degenerate cone at each vertex of Bh for all sufficiently 
small non-zero e. 

3. The condition: ¥= in (|4.1[) is independent of the normalization of the 
eigenf unction, $i. 

Proof of Theorem By Symmetry Remark 12.11 we may without loss of generality 
consider the specific Bh vertex: K, = K. The local character of all others is identical. 

We consider a perturbation of K, K + k, with |k| small. We express $ G 
as $(x;k) = e''''^-0(x; k), where '0(x;k) is A- periodic. The eigenvalue problem for 
k = K + K takes the form: 

( -(Vx + ^(K+«;))^ + y(x) ) V'(x; K + At) = /i(K + k)V'(x; K + k) , (4.4) 
V'(x + v; K + k) = ■0(x; K + k), for aU v e A . (4.5) 

Let /io = M^"-* = m(K) be the double eigenvalue and let V'*'°-' be in the correspond- 
ing two-dimensional eigenspace. Express /i(K + k) and ?/'(x; K + k) as: 

^(K + k) = + V(x;K + k) = + (4.6) 

where V'*'^'' is to be chosen orthogonal to the nullspace of H(K.) — fi^^^ I, and ^j,^-^^ are 
corrections to be determined. Substituting (|4.6p into the eigenvalue problem (|4.4p - 
(|4.5p we obtain: 

(7J(K) -A^o )^^'^ 

= 2iK • (V + iK) - K • K + ^(1) j ij'-^'> 
+ (2iK ■ {V + iK) - n- K, + ^(1) j V'^"^ 

^ F('\ (4.7) 
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Since ijj^^^ is in the Lp^^. ^- nuUspace of H(K.) — figl, we write it as 

= a0i(x) + (x), where (4.8) 
0,(x) = e-^-$,(x). J = 1,2 (4.9) 

Here (f>i and (/)2 are normahzed eigenstates with Fourier expansions as in part 3 of 
Proposition and a,/? are constants to be determined. 

We now turn to the construction of -0^^^ Introduce the orthogonal projections: 
Q||, onto the two-dimensional kernel of i?(K) — iiqI, and Q± = I — Q\\. Note that 



QllV;(i) = QxV^°) = 0, and Qx^A^^^ = V'^^^ . (4.10) 

We next seek a solution to (|4.7p by solving the following system for ip'^^^ and /x'^-'^): 

(i?(K)-A^o )V'<'^ = gx i^(^)(a,/3,^,A^(^\V^^^) (4.11) 
= g||i^(i)(a, /.(I), V^i)) (4.12) 



Equation (|4.12l) is a system of two equations obtained by setting the projections 
of F'^^^ onto (j)i and 02 equal to zero. Our strategy is to solve (|4.1ip for ■0*^^^ as a 
continuous functional of a, /3, k, /i^^-* with appropriate estimates, then substitute the 
result into (|4.12p to obtain a closed bifurcation equation. This is a linear homoge- 
neous system of the form A1(/i*^^\ K)(a, /?)* = 0. The function k ^ ij,^^\k) is then 
determined by the condition that det 7W(/i^"'^) , k) = 0. 

Written out in detail, the system (|4. lip - (14. 121) becomes: 

(iJ(K) -1X0)^(^1= Qx (^2iK-{V + iK) - k • k + ^^^^^ V^^^ 

+ Qx (2iK • (V + «K)) (4.13) 

Ql\ (^2iK- (V + iK) - K + 0(°) 

+ Q|I ( 2iK- (V + iK) )V^i) = (4.14) 

Introduce the resolvent operator: 

R-Kif^o) = ( H{K) - fio ly' 

defined as a bounded map from Q±L'^ to Q±H^(M.'^/Ah). Equation (|4.13p for ijj^^^ can 
be rewritten as: 

(/ + R-k{hq)Q± ( -2iK-(V + iK) + K-K - M^^^ ) ) 0^^^ 
= R-Kif^o) Q± i^iK ■ (V + zK)) (4.15) 

In several equations above we have used (|4.10p . 
By elliptic regularity, the mapping 

f^Af^ i?K(Mo)Q± ( -2iK • (V + zK) + K • K - Ai^i) ) / 
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is a bounded operator on i/'*(M^/A;,), for any s. Furthermore, for + | sufficiently 
small, the operator norm of A is less than one, (/ + A)^^ exists, and hence (I4.15|) is 
uniquely solvable in Q±H'^{R'^ /Ah): 

V'^^^ = + RkMQ± ( -2m-(V + iK) + K-K - //(I) 
o i?K(/^o) Q± (2iK- (V + zK) 

Since is given by (|4.8p . ip^^'^ is clearly linear in a and /? and we write: 

^(1) = c(i)[«;,Ai(i)](x) a + c(2)[k,^(i)](x) /?, (4.16) 

where (K,/i^^-*) i~> c'--'-*[k, /i^^-*] is a smooth mapping from a neighborhood of (0,0) e 
R2 X C into iJ2(M2/A/,) satisfying the bound: 

||c(^)1|h^ *SC(|«:| + |/i«|), j = l,2 . 

Note that Q||c(^) = 0, j = 1,2. 

We next substitute (|4.16p into (|4.14p to obtain a system of two homogeneous 
linear equations for a and /3. Using the relations: 

Vk<^, = = e-*K-^V$„ <0„ 0,> = <$„ $,> = , j = 1, 2 

C(-'')KAi(i)](x) =def e^^-V^)K/.(i)](x), = 0, i,j = l,2 (4.17) 

we have: 

Mifi'-^ln) i 1 = 0, (4.18) 
where M.{^^^\ n) is the 2x2 matrix given by: 



2) 



+ 



^t(l) - K ■ K + <'I>l,2iK ■ V<I>1> <'I>1,2«K • 

< <J>2, 2jK ■ V^l > - K ■ K + < 'I>2, 2«K ■ V'I>2 > 

< $i,2m- VC(i){k,m(i)) > < $i,2m- VC(2)(K,/i(l)) > 

< <I>2,2iK- VC(i){K,At'^)) > < <I'2,2jK- VC(2)(K,/i(i)) > 



Thus, /i(K + k) = /x'-"'' + ^S'^^ is an eigenvalue for the spectral problem (|4.4p - (|4.5p if 
and only if /i^^^^ solves: 

detX(/i(i),K) = 0. (4.19) 

Equation ()4.19p is an equation for /i*-^-*, which characterizes the splitting of the double 
eigenvalue at k 7^ 0. We now proceed to show that if the nondegeneracy condition 
(|4.ip holds, then the solution set of (|4.19p is locally conic. 
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We anticipate that a solution n^-'-^ = 0(1^1) and hence C^^^ = 0{\k\). This 
motivates expanding M. as: 

M{fi^^\K)=Mo{fi'^^\K) + Mi{fi^^\K), where (4.20) 
2^<$„..VcI>,> ^« + 2z<ci>„..V<i>.> j' 

(4.21) 

A^M,(/i«,A^) = 0(|«:|-|M^'^l + \^\') (4.22) 

Now we claim that, in a neighborhood of k = 0, the solutions of (|4.19p are well 
approximated by those of the truncated equation: 

detMoii^.n) = det ^. ' \ n- /I vri \ =0. (4.23) 

^ ' ^ \^ 2i<$2, K ■ V<i>i> 1/ + 2z<$2, • V<i>2> / ^ ^ 

We shall first solve (I4.23P fProposition l4.2l) and then show that the solutions of (|4.19p 
are small corrections to these (Proposition I4.3p . 



Proposition 4.2. 

V$a>= 0, a= 1,2 . (4.24) 



2i <$i,K-V$2>= 2i <$2,K-V$i> 



-3 area(17) ( ^ c(m)2 f ] 



X (ki + iK2) 



= -Aj X (ki + ZK2) ; see (|i?T|) . (4.25) 
We prove Proposition 14.21 iust below. A consequence is that AAq simplifies to 

Moi.;.) = ( , \ . ^ (^1 + ^^^) ] (4.26) 

and therefore 

det7Wo(i^;«) = - |Ajp |«:|2 



3 X area(il) x ^ c(m)^ 



me 5 



2 

2 



X (4.27) 



Therefore, the truncation of to Mo, yields det A4o{h', k) = and a locally conical 
dispersion relation, provided Aj 7^ 0. 

Proof. [Proof of Proposition 14. 2| Recall that $1 e ^ and $2 e r ^'"'^ given 
by p.6ip and (12.621) . respectively. We first consider the diagonal elements and claim 
($1, K ■ V$i) = 0. To see this, first apply k • V to $1, given by (|2.63p and obtain: 

= iK - c(m) ( K:"e*^"-^ + f RK^e'"^^ '' + r R^K^e'^'^^ ''^ . 

me 5 
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Therefore, since e~*^"" e'^^^dx = if m 7^ n we have 
K ■ <$i, V$i> = area(ri) ^ \c{m)\^ k ■ (K^ + \t\^ RK^ + \t\^ R^K^ ) = 0. 



me 5 



(4.28) 



The latter equality holds since jrp = 1, / — i? is invertible (spec(i?) = {t,t}) and 
I + R + R^ = = {I-I){I-R)-^ = 0. Similarly, <$2,k- V$2> = 0. 

Thus, we have shown that the diagonal elements vanish, (|4.24p . 



One can check directly that the off-diagonal elements satisfy A^o,i2 = ■M.o.21- 



2i<$i,K- V$2> = 2i<$2,K- V$i> . (4.29) 
Furthermore, using (|2.6ip . (I2.62p we have 

2i<$i,K- V-I>2> 

= -2 X area(r!) ^ ^{m}^ k ■ (K™ + RK^ + f ' R^K^ ) , 



= -2 X aveain) ^ c(m)2 k • + r i? + (r i?)^ j . (4.30) 

Note, by (I2.28P that tR has an eigenvalue f with corresponding eigenvector ( = 
2~^/2(l,i)* and an eigenvalue 1 with corresponding eigenvector ( = 2~^/2(l, — i)*. We 
express K™ as 

= <C,K™>C + <C,K->C, 
where for a, b e we defined (a, b) = a • b. Clearly, 

(/ + T R+{t i?)^) C = 3 C and (/ +t R+{t i?)^) ( = (l + f + (f)^) C = 0. 

Therefore, 

K • (/ +T R+{t Rf^ K™ = 3 <C, K™> K . C = 3 (C • K™) x {k ■ () (4.31) 
Substitution of (j43T|) into (|4:30| we obtain: 
2i <$i,K- V$2> 



-6 X area(r2) ^ c(m)2 (C • K™) x (k • C) 
me 5 



—3 X area(r2) ^ c(m)2 



me 5 



1 - .K^ 



1 

K2 / V -« 



me 5 



3 X area(f^) ^ c(m)2 ( M.K'" X {KI+IK2), (4.32) 



which by the definition of Aj in (|4.ip implies (|4.25p . Thus, detA^o(z^,K) = i^2_|;^jj|2|^|2 
which proves (|4.27l) . Furthermore, the solutions of detA^o('^iK) = define a non- 
trivial conical surface provided: 

Aj = 3 X area(f7) x ^ c(m)2 ( M • K"" 7^ 0. (4.33) 

me5 ^ * ^ 
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The proof of Proposition is complete. □ 

To complete the proof of Theorem 14.11 we next show that the local character of 
solutions to (|4.19p is, for |k| small, essentially that derived in Proposition 14.21 for the 
solutions of (1423)) . 

Proposition 4.3. Suppose Aj, defined in (|4.33p . is non-zero. Then, in a neigh- 
borhood of any K or K' point, the dispersion surface is conic. Specifically, the eigen- 
value equation det A4{ii^^\ n) = (see (14.191) ) defines, in a neighborhood U cz R"^ of 
K = 0, two functions: 

= |A„| \k\ {1 + E+{^)), ^i'•l\>,) = -|A„| |«| {l + E^{^)), (4.34) 

where E+(k) — > as k — > and E+(k) is Lipschitz continuous in k. 
Proof By (|4.19|) . and Proposition 14.21 satisfies: 

(m^i))' = \Xtf l^f + .921 +512 +503(«;), (4.35) 

where grs are smooth functions satisfying the bounds: 

I grsiti,^^) I < C l/ir I'sl' 

for ^ 1, |k| ^ 1. We now construct ^^_^\n). The construction of ^'"^ (k) is similar. 

Set ^^^^ = |Ajj| |k| (1 + 77). Substitution into (|4.35p and using that Aj 7^ 0, we find 
that 77 satisfies: 

G{r],K) = 2rj + rj^ + Ji{rj,K) + J2{n) = Q. 

Here Ji and J2 are smooth functions of r; and Lipschitz continuous functions of 
K, such that: Ji{r],K.) = 0{\k\), d,^Ji{T],K) = 0(|k|), J2(k) = 0{\k\) as |k| 0. 
Thus, G(r], k) and drjGij], k) are Lipschitz continuous in (77, k) with G(0, 0) = and 
driG{0, 0) = 2 7^ 0. It follows easily that there exists rj = E{k) defined and Lip- 
schitz continuous in a neig hborhood C/ c of k = 0, such that £:(0) = and 
G{E{k),k) = for all keU.Q 

5. Main Theorem: Conical singularity in dispersion surfaces. Assume 
that y is a honeycomb lattice potential in the sense of Definition 12.11 Since V e 
C"^ {M."^ / Afi) , its Fourier coefficients satisfy 

Vel\Z^), I.e. II y = \Vrr.\<^. (5.1) 

meZ2 

Theorem 5.1. Conical singularities and the dispersion surfaces of H^"^^ 
Let V{'x) honeycomb lattice potential. Assume further that the Fourier coefficient 
of V , Vi^i, is non-vanishing, i.e. 

Vi,i = [ e-*('=^i+'=^)y V{y) dy ^0 . (5.2) 
Jn 

There exists a countable and closed set C c M such that for any vertex K* of Bh and 
all e ^ C the following holds: 
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1. There exists a Floquet-Bloch eigenpair $'^(x; K,), /i'(K») such that 

yLt'^(K,) is an ^ - eigenvalue of H^^^ of multiplicity one, with corre- 
sponding eigenf unction, $'^(x;K*). 

/i'^(K*) is an ^ - eigenvalue of H^^^ of multiplicity one, with corre- 
sponding eigenf unction, $'^(— x;K»). 

/i'^(K») is not an ^- eigenvalue of H^^\ 

2. There exist 6^ > 0, > and Floquet-Bloch eigenpairs: (<i>'^(x; k), /i^(k)) 
and ($l(x; k), /xl(k)), and Lipschitz continuous functions, £'+(k), defined 
for |k— K*| < Se, such that 

/x;(k)-/x-(K,) = + |A^| |k-K.| (1 + Elik) ) and 

f,%{k) ^ t,%K.) = -\XI\ |k-K.| (1 + EL{k)), 

where ¥= is given in terms of "I>'^(x;K*) by the expression in (|4.ip 
and |-E+(k)| ^ C^jk — Kj,|. Thus, in a neighborhood of the point (k, ^) = 
(K»,/x^) G M.^, the dispersion surface is conic. 

3. There exists £° > 0, such that for all e e (-e°, £'')\{0} 

(i) eVi^i > =^ conical intersection o/ I'** and 2"'* dispersion surfaces 

(ii) eVi^i < =^ conical intersection of 2^'^ and 3^"^ dispersion surfaces . 

Remark 5.2. Part 3 of Theorem \5.1\ gives conditions for intersections of the first 
and second band dispersion surfaces or interesections of the second and third. As the 
magnitude of e is increased it is possible that there are crossings among the ^ - 
eigenvalues of H^^\ so in general the theorem does not specify which band dispersion 
■surfaces intersect. 

5.1. Outline of the proof of Theorem 15.11 By Symmetry Remark 12.11 it 
suffices to prove Theorem 15.11 for K« = K. We have seen that the central pomt is 
to verify for all e, except possibly those in a closed countable exceptional set, that 
hypotheses (hl-h3) of Theorem 14.11 hold. These hypotheses state that H^'^^ has simple 
^K, r ^^'^ ^K, T eigenvalues which are related by symmetry, which are not 
eigenvalues, and moreover that A| ^ 0. We proceed as follows. 

In section [5] we show that there is a positive number, £°, such that for all e e 
(-e°,£°)\{0} (hl-h3) of Theorem O hold. That is, the conclusions of Theorem O 
hold for all sufficiently small, non-zero e. In section FTSl we introduce the key tool, 
a renormalized determinant, to detect and track the ^ eigenvalues of i/^*^-* for 
cr = 1,T, f. A continuation argument is then implemented using tools from complex 
function theory in section[51 to pass to large e. We now embark on the detailed proofs. 

6. Proof of Main Theorem 15.11 for small e. We begin the proof of Theorem 
15.11 by first establishing it for some interval — e° < e < e^, where is positive but 
possibly small. We shall consider the eigenvalue problem for iJ^^^ on the three eigen- 
spaces of TZ: iK.,T'-^K,,r ^^^^ Lj^^y. 

i/^-^) $(x;K.) = [-A + eV{x)] $(x;K*)= ^(K.)$(x;K.) (6.1) 
4>(x + v,K*) = e*^*-^$(x,K«), x6R2 
7e[$(-;K.)] (x) = a$(x;K,), where CTe{l,T,f}. 

An eigenstate <&(x; K,) in ^ is, by Proposition l2.61 of the form: 
$(x;K,) = c(m;$) ( e*^" " + ct e''^^" " + a e^^'^" '') . (6.2) 

me5 
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The summation is over the set, S, introduced in Definition l2.4l Note that by Proposi- 
tion l2. Gl and Remark [2.3[ solutions to the eigenvalue problem on ^ can be obtained 
from those in ^ via the symmetry: $(x) i-^ $(— x). 

Recall that c$(m) or c(m; $) denote the ^- Fourier coefficients of <&. Our 
next task is to reformulate the eigenvalue problem (|6.ip as an equivalent algebraic 
problem for the Fourier coefficients {c(m; $(•; Kt))}ine5- First, applying —A — /i to 
$, given by (j6.2|) . and using that R is orthogonal, we have that 

(-A-^)<f(x;K.) 

me 5 

(6.3) 

Next, we claim that V^(x)$(x; K«) e ^. Indeed, since is TZ- invariant, 7?,[V^](x) = 
V{R*x) = V{x). Moreover, since $(-;K*) e ^, we have ??.[$] = cr$. Therefore 

7^[^/$] = T/(i?*x) $(i?*x;K,) = V^(x) a $(x;KO = a 

Therefore, by Proposition 12.61 l/^$(-;K*) has the expansion 

y(x)$(x;K,) = c(m;F$) ( e*^" '' + e'-"^"'^ + e^-"'^"'^) , (6.4) 

c(m;F$) = p ^ e-'K^ y ^(y)$(y;KO dy . (6.5) 

Furthermore, with the notation qk • x = (^iki + (72^2) • x, 



c(m; F$) 



■^V<f){y) dy 



X j dy |^e*(K^.-K?'+qk) y _^ ^ gi(J?,K;-K^+qk) y ^ ^ gj(i?,2K:-K^+qk) y 



X l^e'^'l"^™"''))'')'^ + a p'(q-(m-'^r))k)-y _^ ^ g»(q-(m-TC"r))k)-y" 

E yqc(r;$) 

X [5(q-(m-r)) + 5 (q - (m - 7^r)) + ct (5 (q - (m - 7^2r))) ] 
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Thus, 

c(m;F<f) = 2 /C,(m,r) c(r;<i>), 

where (recall (I2.5ip ') 

/C^(m,r) = Vm-r + ^ Vm-7Zr + ^ V^-TC^r (6-6) 
^'^1 — ri ,7n2 — r2 ^" '-^ ^^ii +r2 ,m2 +r2 — ri — 1 ^" '-^ ^^ii+ri — r2 + l,rM2+ri ■ 

Summarizing, we have 

Proposition 6.1. Let a e {1,t, f}. T/ien, the spectral problem (16. ip on L"^^ ^ 
is equivalent to algebraic eigenvalue problem for c(m) = c(m; $) anrf /i; 

(|K™|^-^) c(m) + £ XI ^'T(m,r)c(r) = 0, meS, (6.7) 

w/iere {c(m)}me5 e ^'^(5). 

To fix ideas, let K« = K; starting with K', we would proceed similarly. For e = 0, 
we have the algebraic eigenvalue problem: 

+ miki + m2k2|^ — c(m) = 0, meS. (6.8) 

Equation (|6.8p . viewed as an eigenvalue problem for ({c(m)}^gz2, A*) is equivalent 
to the eigenvalue problem for —A on treated in Proposition 13.11 Restated in 
terms of Fourier coefficients. Proposition 13.11 states that /i^"^ = |Kp is an eigenvalue 
of multiplicity three with corresponding eigenvectors: 

|K|^ 

|K|2 = |i?K|2 = |K + k2|2 

|K|2 = |i?2K|2 = |K-ki|2 

Recall from Definition \2.4\ that the equivalence class of indices {(0, 0), (0, 1), (— 1, 0)} 
has as its representative in S the point (0, 1). 

The eigenvalue problem ()6.8p has a one dimensional ^- eigenspace with eigen- 
pair: 

= |K + k2p = |Kp, c(mi,TO2) = (Smi,m2-i, (mi,r7i2)e5 
corresponding to the ^ eigenstate of Hq : 

= e'(K+k2)-x ^ -g«(K-ki)-x ^ ^g.K.x 

= CT e'^ '' ( 1 + a e*'^^ '' + ae"'''! '') 

We seek a solution of (|6.7p for e varying in a small open interval about e = 0. We 
proceed via a Lyapunov-Schmidt reduction argument. First, decompose the system 
(|6.7p into coupled equations for: 

C||=c(0,l) e C, and {cj_(m)}^e5^ e /^(5-^) , (6.9) 



C^i°^(mi,m2) = Smi,m2 
4°'' (mi, 777-2) = <5mi,m2-l 
4°)(-l,-2)=5™.+l,™2 



,(0) 



,(0) 
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where 



and rewrite (j6.7p as a coupled system for C|| and c±: 

[{K^'^f-fi + £ /C,(0, 1,0,1)] c|| + e 2 ^a(0,l,r) cx(r) 



(6.10) 



0, (6.11) 



£ /C,(m,0,l) C|| + |K"f -M cx(m) 



+ £ 2 /Ccr(m,r) cx(r) = 0, m e iS"" 

re5^ 



(6.12) 



We next seek a solution of (|6.1ip - (|6.12p . for e small, in a neighborhood of the solution 
to the £ = problem: c?. = l,fi^°'> = |K|2,cx(r) = 0,r e 5-^. 



We begin by solving the second equation in (|6.12p for c± as a function of the 
scalar parameter C|| . For e small, the operator to be inverted is diagonally dominant 
with diagonal elements: |K™|^ — /i, which we bound from below for m e 5^. By the 
relations (|2.22p we have 



|K'= 



|K + TOiki + m2k2| — fi 
|K|^ — ^ + {m\ + 1712 — "117712 + TOi — TO2) , q 



47r 
a^/3 



If fj, varies near ^^^^ = |K|^, then 



IK" 



5= me S 



for some ci > 0. We now rewrite the equation for c± as: 



Sm.r + 



2 /C^(m,r) 



cx(r) 



- £ cii 



/CCT(m, 0, 1) 



|K>"|2 

= £ C|| F^{n), meS^ 

or, more compactly, 

(/ + sTKAf^) )cx = £C|| F^(Ai) 
Recall Young's inequality, which states that the operator defined by 

nfim) = 2i(m,r)/(r) 



(6.13) 
(6.14) 



satisfies the bound 



\\TLf\\p{s^)) < Cl where 

Cl = sup 2 |i(m,r)| + sup2|L(m,r) 



(6.15) 
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We apply (|6.15p with L(m,r) = !C„{in,r), defined by (|6.6p . and conclude using the 
bound 



||K™|2-/i| ^ l + |m|2^ 
and that V = {Vin}mes^ ^ ^^('5^) (recall (|5.1I) ). that the operator 

maps /2(5-'-) '2('^^) ■^ith the bound 

II rK.(M)/ ILi(s^) =S C ||/||p(z2). (6.16) 

Here, = I]™e5-(1 + |mn2|/(m)p. 

Proposition 6.2. There exists > such that for all \e\ < e" and any f e 

{I + eTkAi^)) = { , (6.17) 

has a unique solution cx = e I'^iS'^), analytic in e, satisfying \\c±\\i2(^g±^ ^ 
2||f|l;|(5i) ■ 

We now apply Proposition l6.2l to solve (|6.14p to obtain 



(r) = £C|| [{I + eTKA^^))-'F'^{^^)]{r). (6.18) 



Substitution into (16. lip yields a closed scalar equation for C|| of the form Ad^in, e)c[| = 
0, which has a non-trivial solution if and only if: 

EE \Kf-^i + e/C40,l,0,l) + 2 X;40,l,r) [(/ + eTk^ (m))"' F"(^)] (r) = 

(6.19) 

Ma{fi,e) is analytic in a neig hborhood of (/i,£) = (Ai(°\0) = (|K|2,0). Clearly, 
Mcr{fi^°\0) = and (5^M(j(^°,0) = —1. By the implicit function theorem, there 
exists > such that defined in a complex neighborhood of the interval |e| < e°, 
there is an analytic function e ^ ^'^ , such that 

Ma{^J',e) = 0, for - e° < £ < e° 

Thus, we take C|| = 1 and via (|6.18p - (|6.19p our solution for |e| < e*^ is 

^= = |K|2 + e/C^(0,l,0,l) + C'(e2) 

C\ = C(0, 1) EE 1 



where F^,M = -^^T^T^' ™ ^ '^^ ' (6-20) 
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From the definition of /Ccr(m,r), displayed in (|6.6p . we find: 
/C,(0,l,0,l) = Fo,o + al^i,! + aFo,i = Fo,o + Fi,i(a + a), a=l,r,f. (6.21) 
The latter equality uses: 

(a) constraints on Krei,m2 by TZ— symmetry of V (Vo._i = Vi.i) and that 

(b) V is even (Vb,-i = Vb,i). 

Furthermore, Vi_i is real, since T^(x) is even and real {Vi^i = V-i-i = Vi^i). There- 
fore, 1C„{Q, 1, 0, 1) is real, as expected. 

The small e perturbation theory of the three-dimensional eigenspace is now sum- 
marized: 

Proposition 6.3. Assume 0. Then, there exists e° > such that for 

< |e| < the multiplicity three eigenvalue n = |Kp perturbs to 2- dimensional 
and l-dimensional eigenspaces with corresponding eigenvalues /i^(K) and /i'^(K) as 
follows: 

1. /i^(K) is of geometric multiplicity 2 with a 2-dimensional eigenspace X,- © 
Xt= c ^ © ^ (^iwen fey; 

fi'{K) = |K|2 + £ (l/o,o + 2Fi,i cos(27r/3)) + 0{e^) 

= |K|2 + £(yo,o-V^i,i) + 0(£') (6.22) 

mi/i eigenstates $^ e -/^^ ^ and <i>| e ^kt' obtained by the symmetry (see 
RemarklEM): 

$2(x;K) = $i(-x;K), 

wii/i Fourier expansions: 
$f(x,K) = 2 c"(m) (e'^" " + f e*«^'" " + r e'^'^'" ") . (6.23) 

me 5 

$§(x,K) = 2 (e''^""' + r e*-^^"-^ + f e*'"'^'" '') , and 

me5 

(6.24) 

^. /i'^(K) zs a simple eigenvalue with eigenspace Xi c ^: 

M^(K)= |K|2 + e(Vb,o + 21/i,i) +0(£2) (6.25) 
with eigenstate $; 

l'^(x,K) = 2 5'^(m) (e*^"-^ + e'-"^"'^ + e^^'^'"'') . (6.26) 

me 5 

Proposition 16.31 implies that the double-eigenvalue hypotheses of Theorem 14.11 
holds for e positive and small. In particular, by (|6.22p and (|6.25l) 

If eVi,i > 0,then 

(K) = ni'^ (K) < (K) < mI') (K) . . . (6.27) 

and if eVi.i < 0, then 

(K) < (K) = (K) < mI') (K) < . . . . (6.28) 
By Theorem 14. 1[ assuming ^ 0: 
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(i) if eVi i < 0, the dispersion surfaces k i-^- /xi(k) and and k i-^- /i2(k) intersect 
conically at the vertices of Bh 

(ii) if eVi^i > 0, the dispersion surfaces k i— > /^2(k) and and k i— > /i3(k) intersect 
conicahy at the vertices of Bh ■ 

So, in order to apply Theorem 14.11 it remains to check that A| 7^ . Here, A| is the 
expression given in (|4.ip . For e small we have 

Aj = 3 area(r2) 

where we have used that 0*^(0, 1) = 1, (|6.20p and that ||c^||i2(5X) = 0{e). Note also 
that 

Therefore, for \e\ < e^, with chosen sufficiently small, 

2 

lA^I^ = 16 area(l])2 ^ + 0(e). (6.30) 

This completes the proof of our main theorem. Theorem 15.11 for the case where e is 
taken to be sufhciently small. We now turn to extending Theorem 15. II to large e. 

7. Characterization of eigenvalues of H^'^^ for large e. To extend the as- 
sertions of Theorem 15.11 to large values of e, we introduce a characterization of the 
^- eigenvalues of the eigenvalue problem (16. ip as zeros of an analytic function of 

e. 

Since we can add an arbitrary constant to the potential, by redefinition of the 
eigenvalue parameter, /x, we may assume without loss of generality that 

V{k) s; Knax. 

Assume first that e e C and > 0. Then, H'^"'^ - fil = -A + eV - fil = 
(—A + eV + /) — (/i + 1)/. The eigenvalue problem (|6.1I) may be rewritten as 

(-A + £y + /)$-(/i+l)$ = 0, ueL^,^. (7.1) 

Now for any real e > we have — A + eV + 1^1. Hence we introduced 

r(e) = {I-A + sV)-\ (7.2) 

which exists as a bounded operator from ^ to H-^ ^ and obtain the following 
Lippmann - Schwinger equation, equivalent to the eigenvalue problem (j6.ip : 

[/ - ifi+l) T{e) ] $ = 0, $eL^,, . (7.3) 

We now show that if 5Re < 0, we also obtain an equation of the same type as in 
(|7.3p . In this case, we observe that e {V — T4iax) 5= 0. Therefore, —A + e{V — Knax) + 
I ^ I and we rewrite (|6.ip as 

(-A + eiV- K.ax) + + eKnax) * = 0, $ 6 L^^,, . (7.4) 

t For !R£ > and / smooth, we have Jf {(— A + eV + /) /, /) > ||/P- Hence the nuUspace of 
—A + eV + / and its adjoint are {0}. By elliptic regularity theory — A + eV + / is invertible on ^. 



+ 0{e), 



(6.29) 
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If for e < we define T(e) = 1— A + e{V — Knax), then (j6.ip is equivalent to 

[ / - + 1 - eKnax) f{e) ] $ = 0, $ e L^^^ . (7.5) 

For the remainder of this section we shall assume > and work with the form 
of the eigenvalue problem given in (|7.3p . The analysis below applies with only trivial 
modifications to the case e < and the form of the eigenvalue problem given in (|7.5p . 

For each e > 0, we would like to characterize ^- eigenvalues, /i(£), as points 
where the determinant of the operator / — (^ + 1) T(e) vanishes. To define the deter- 
minant of / — zT, one requires that T be trace class. Although T{e) is compact on 
^K.cn it is not trace class. Indeed, in spatial dimension two, Xj, the j*'' eigenvalue of 
~Ak + W acting in L^^^j^ satisfies the asymptotics Xj ~ \j\ (Weyl). Therefore 

trace(T(£)) = Y,\X,\-' ~ = O) . 

3 J 

The divergence of the determinant can be removed if we work with the regularized 
or renormalized determinant; see [3 [13]. Note that T{e) is Hilbert Schmidt, i.e. 

3 3 

For a Hilbert-Schmidt operator. A, i.e. tr(yl^) < cx), define 

R^{A) ^ [I + A]e-^-I. (7.6) 
Note that I + A\s singular if and only if / + R2{A) = {I + A)e^^ is singular. 

Since = 1 — z — (s — l)e~*^(is, we have 
R^{z) = (1 + z)e-^ - 1 = -z2 (^1 + (1 + z) /q (s - l)e-'^ds^ . Therefore 

i?2(A) = -A^ A) j\s - \)e-'^ds^ . (7.7) 

Since A^ , is trace class and the second factor is bounded, R2{A) is trace class. There- 
fore the regularized determinant oi I + A: 

det2(/+A) = det(/+i?2(A) ), (7.8) 

is well-defined. With A= -{^,+ l)r(e), we have the following [Tl iniITT] : 
Theorem 7.1. Let a take on the values 1,t or f. 

1. £ t—^ T{e) is an analytic mapping from {e e : Ke > 0} to the space of 
Hilbert-Schmidt operators on ^. 

2. ForT^e), considered as a mapping on L^^, define: 

f,(M,£) = det2(/-(M+l)T(£)). (7.9) 

The mapping (/i,e) i— >■ Saif-j^), which takes (/i,£) e (IRe > 0) to C is 
analytic. 
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3. For £ real, jjL is an ^- eigenvalue of the eigenvalue problem (|6.ip if and 
only if 

£Al^,e) = 0. (7.10) 

4- For e real, fi is an ^ eigenvalue of (|6.ip of geometric multiplicity m if 
and only if fi is a root of (j7.10p of multiplicity m. 

8. Continuation past a critical e. In section [S] we proved Theorem 15.11 for 
all e e (— e", e°)\{0}, with > sufficiently small, by establishing the following 
properties: 

I. /Lt^(K») is a simple ^ eigenvalue of H^^"^ with corresponding 1 - dimen- 
sional eigenspace Xr = span{ (x; K*) } c ^k^.t- 
II. /i'^(K*) is a simple ^ eigenvalue of H^^^ with corresponding 1 - dimen- 
sional eigenspace Xt= = span{ (^^jK*) } c ^. 

III. /i^(K*) is not a L^^^ eigenvalue of H^^K 

IV. We have 

= 2 c(m,/.^e)2 (])-KT^O, (8.1) 

where c(m, /i^, e) are Fourier coefficients of $i [/x'^(K«), e] (x), an ^ eigen- 
function of i/'-^^ with eigenvalue /j,'^ = /i^(K«) ; see Proposition 12.61 

We next study the persistence of properties I. -IV. for e of arbitrary size. 

8.1. Continuation strategy. Denote by A, the set of all e > for which at least 
one of the properties I. -IV. fail. With e'^ given as above, we clearly have Acz [e", cx)). 
The main result of this section is that 

A is contained in a countable closed set. (8.2) 

Once (|8.2I) is shown, we'll have completed the proof of Theorem 15. 11 our main result. 

Our continuation strategy is based on the following general 
Lemma 8.1. 

Let A c: (e'^,cx)) with e*^ > 0. Then one of the following assertions holds: 

(1) A is contained in a closed countable set. 

(2) There exists Ec e (0,oo) for which the set A n [0,£c) is contained in a closed 
countable set, but for any e' > £c, the set A n [0,£') is not contained in a closed 
countable set. 

The main work of this section is to prove (|8.2p for ^ = ^ by precluding option 
(2) of Lemma [8. II This suggests introducing the notion of a critical value of e: 

Definition 8.2 (Critical £c). Call a real and positive number £c critical if there 
is an increasing sequence {e^} tending to £c and a corresponding sequence of geometric 
multiplicity-two L^- eigenvalues, {piv}, such that 

(a) properties I. -IV. above, with £ replaced by and /i^ replaced by fi^, hold for 
all V = 1,2, ... , and 

(b ) for £ = £c and pic = m"^^ = lim^_>oo fJ.^ < X> at least one of the properties 
L-IV. does not hold. 
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To prove Lemma l8.ll we use the following: 

Lemma 8.3. Let = ei < £2 < ■ • ■ , and let eo-j = Imii^^oc £i^- (Perhaps Src = co.) 
Suppose A n [0,Si,) is contained in a closed countable set C^, Jar each > 1. Then, 
A n [0,£oo) is contained in a closed countable set C. 

First let's use Lemma 15751 to prove Lemma [01 We then give the proof of Lemma [5751 
Proof of Lemma \8.1i 

Let £c = sup{£ e (0, 00) -.An [0,e) is contained in a closed countable set.}. Clearly 
< £" £c 00. If £c = 00, then option (1) holds, thanks to Lemma 18.31 And 
if £c < 03, then by definition, A n [0,£') is not contained in a closed countable set 
for any e' > £c- Again applying Lemma [8.31 shows that A n [0,£c) is contained in a 
countable closed set. In this case, (2) holds and the proof of Lemma [8. II is complete. 

Proof of Lemma \8.S\ Define 

C=[j{C^r^ [£^_i,£^] ) u {£^ : 1/ > 0} 

and set C = C if £00 = 00 and C = C I if £00 < CO. One checks easily that 

A D [0, £00) c C, C is countable, and C is closed. This completes the proof of Lemma 

[Ql 

Remark 8.1. Our proof will in fact show that defined on is a function 
£ 1-^ /3(e) J analytic on the connected components of the complement of a countable 
closed set C and such that: (a) I3{e) is a multiplicity two L^^- eigenvalue of H'^^^ and 
(b) A| = A|j(/3(£),£) 7^ 0, outside C. 

In section 15721 we develop general complex function theoretic tools concerning the 
variation of zeros. A, of a real-analytic function, P(A, z) as a parameter, z, is varied. 
This result is later applied to P = £^{fi,£) for a = 1 and tr = r, the renormalized 
determinant whose zeros are ^ eigenvalues of iJ^^^; see section[7l In section [873l we 
study the general problem of constructing non-trivial null- vectors oi NxN matrices of 
rank A^— 1, which depend analytically on the matrix entries. The results of sections [8. 2 1 
and l8.3l are applied to the construction of an eigenfunction of H^^^ and a corresponding 
simple ^- eigenvalue, /?(£), which vary analytically with e. It follows via (|8.ip that 
A| = X^{f3{e), e), varies analytically with e. These results play a central role in the 
continuation strategy outlined above. 

8.2. Picking a branch. Let 

U = {(A, z) 6 C2 : |A| < £1, \z\ < £2} (8.3) 

where £1 and £2 are given positive numbers. Suppose we are given an analytic func- 
tion P : U ^ C and an analytic mapping F : U ^ C™. We make the following 
Assumptions: 

(Al) If (A, z) e U, P(A, z) = and z 6 M, then A e M. 

(A2) There exists {(A^, z^)} c [/, 1/ ^ 1 tending to (0, 0) as ^ 00, such that for 
each iy^l,z^e M\{0}, P(A^,z^) = 0, ^aP(A^,z^) ^ 0, F(A^,z^) ^ 0. 
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Remark 8.2. With the above setup, we have centered the analysis about (z, A) = 
(0,0). We shall apply the results of this section to an appropriate analytic function 
of{fi,e) centered about {^c,£c)- 

Under assumptions (Al) and (A2) we will prove the following 

Lemma 8.4. There exist 5 > and a real-analytic function l3(z), defined for 
z e (0, 6), such that for all but at most countably many z e (0, S) we have: 

PW{z),z) = 0, dxP{(3{z),z) ^ 0, F{^{z), z)^0 . (8.4) 

Moreover, \imz^Q+ P{z) = 0. 



Proof of Lemma \E^\ Assumption (Al) implies that A i-^ P{\, 0) is not identically zero. 
By the Weierstrass Preparation Theorem !9i , we may write P{\^ z) = H{\, z) ■ P{X, z) 
for (A, z) in some polydisc U = {|A| < £3, \z\ < £4}, where P is a Weierstrass 
polynomial (see (|8.5p below) and iJ is a non- vanishing analytic function. Assumptions 
(Al), (A2) hold also for P, F, U . Moreover, the conclusion of Lemma [8.41 for P, F, U 
implies the conclusion for P,F,U. Therefore, it is enough to prove Lemma 15^ under 
the additional assumption that P is a Weierstrass polynomial. Henceforth we make 
this assumption. Thus, we have for some D ^ 1: 

D-l D 

P(A, z) = A^ + 2 9A^)^' = n - "''(^)) ' (8.5) 

j = i/=l 

where ai{z), . . . , a£)(z) denote the roots of A >— >■ P(A. z) (multiplicity counted), where 
aj(0) = linr aj(z) = 0, j = l,...,D. (8.6) 

Moreover, gj{z) are analytic in \z\ < 64. Note that D ^ 1, since Assumption (A2) 
implies P(0,0) = 0. For k ^ 1, define 

[ £>, for fc = 1 

Qk{z) = < T-rfc f f \ ^ 7 o (8-''') 

The right hand side of (|8.7p is a symmetric polynomial in ai{z), . . . , aoiz) and is 
therefore a polynomial in the coefficients gj{z) of P{X,z) [5], which are analytic in z. 
Consequently, each Qk{z) is an analytic function of z. Moreover, when z is real, the 
av{z) are also real, and therefore, for z e E, Qk{z) 7^ if and only if A i-^ P(A, z) has 
at least k distinct zeros. In particular, for > D + 1, Qk{z) = for all real z, since 
A 1-^ P(A, z) has only D zeros; see (|8.5p . Hence, there exists k with \ ^k ^ D such 
that 

Qk{z) is not identically zero, but Qk{z) = for all k > k. 

Since Qk{z) is analytic on a neighborhood of and not identically zero, there is an 
open interval (0, 61) such that Qk{z) 7^ for all z e (0, i5i). Thus, P(z, A) has at least 
k distinct zeros, for each z e {0,di). On the other hand, Qk^i{z) = and hence 
A 1-^ P{X, z), for real z, never has at least k + 1 distinct zeros. So, A 1-^ P(A, 2) has 
exactly k distinct zeros for z e (0, i5i). We denote these zeros by 



/3i(z) < /32(z) < ••• </3fc(2); 
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they are real by Assumption (Al). Note that each (ik{z) is among the aj{z) {j = 
1, . . . ,D). Hence, by (|8.6p hm2_,Q+ /^^(z) = for each k. 

Fix X e (0, and let mi, . . . , m^. (respectively) be the multiplicities of the zeros 

Pi{x), f32{x), . . . ,f3k{x) of A P(A,z); miH hm^ = D. For z e (0,(5i) close enough 

to X and for each j, there exist rrij zeros of A i— > P(A, z) (multiplicities counted) that 
lie close to fij{x). Unless these nij zeros of A i-^ P(A,z) are all equal, the function 
A I— > P(A, z) would have more than k distinct zeros, which is impossible. Therefore, 
for z e (0,(5i) close to x, and for each j, the polynomial A i-^ P{X,z) has a single 
zero I3j{z) of multiplicity nij, close to Pj{x). In particular, the multiplicities of the 
zeros /3i(z), /32(2;), . . . , /3fe(z) are constant as z e (0,(5i) varies over a small enough 
neighborhood of x. Since x e (0,5i) is arbitrary and since (0,(5i) is connected it 
follows that the multiplicities mi, . . . ,mj, (respectively) of the zeros /3i(z) < /?2(^) < 
• • • < /3fe(z) of A i-» P(A, z) are constant as z varies over the entire interval (0,(5i). 
Therefore, we have 

5 

P(z,A) = Y\{\-p,{z)r^, zG(0,5i), |A|<si (8.8) 
i=i 

Here, fii{z) < P2{z) < ■ ■ ■ < /3j,(z) for each z e (0,(5i) and each nij is a positive 
integer. Now note that each (^j{z) is a real-analytic function on (0,(5i), since /3j{z) is 
a simple zero of A i-^ c/™^ ^P(A, z). 

We now turn to F : ^ C". Let us write F(A, z) = (Fi(A, z), . . . , P™(A, z)). 
For a small positive number p, to be chosen just below, we define: 

G(^) = 7^ S i] fKA,z) ■ JHiyi) ■ ^^'i^l-^^'f^ dxP{\z) dX (8.9) 
27^1 J|A|=p ^(A,2:) 

We can pick p so that P(A, 0) 7^ for |A| = p. Therefore, for small enough 77, if 
|z| < 77, we still have P(A, z) ^ for |A| = p. Fix such p and 77. Then, G{z) is an 
analytic function of z in the disc {|z| < 77}. Moreover a residue calculation shows that 

m 

G{z) =22 Pi (A, ^) • FiCX, z) ■ dxP{\ z) ■ dxPiX, z), (8.10) 

A 1 = 1 

where the sum is over all A in the set: 

{A : |A| <p, P(A,z) = } , 

with multiplicities included in the sum. In particular, if z is real, then the relevant 
A's are also real (see (Al)), and therefore 

m 

G(z) |PKA,z)M5aP(A,z)P (8.11) 

where the sum is over real A e (— p, p) such that P(A, z) = 0. Note that all non-zero 
contributions to the sum ()8.11|) come from A's that are zeros of P with multiplicity 
one. Consequently, for real z, we have G(z) if and only if there exists A e (— p, p) 
such that P(A, z) = 0, dxP{X, z) ^ and P(A, z) ^ 0. 
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Therefore, assumption (A2) teUs us that the analytic function G{z) doesn't vanish 
identically in {|z| < rj}. It follows that we can pick a positive 6, less than mm{ri,di), 
such that 

G(z) 7^ for 0<\z\<S . 

Now suppose z e {0,6). Then, there exists A e {—p,p) such that P{X,z) = 
0, dxP{X, z) ^ 0, F{\, z) ^ 0. This A must be equal to one of the (3j{z), j = I, . . . ,k, 
for which rrij = 1. So, for each z e {0,5) there exists j such that nij = 1 and 
F(/3,(z),z)^0. 

Unfortunately, the above j may depend on z. However, we may simply fix some 
xq 6 (0,(5), and pick jo such that rrijg = 1 and F{l3j„{xo),xo) ¥^ 0. The function 
z 1-^ f3ja{z) is a real-analytic function of z e {0,6). Moreover, we know that the real 
analytic function z >— > F{j3jg{z), z) is not identically zero on (0, 6), since it is nonzero 
for z = xq. So, it can vanish only on a set of discrete points which accumulates at 
or at 6. The proof of Lemma [531 is complete. 

Remark 8.3. We have proven more than asserted in Lemma \8.4\ In fact, 
P{^{z),z) = and dxP{l3{z),z) 7^ for all z e {0,6); and F{P{z),z) ^ for all 
z e (0, 6) except perhaps for countably many z tending to 0. Note that we can arrange 
for this countable sequence not to accumulate at 5 by simply taking 5 to be slightly 
smaller. 

8.3. Linear Algebra. Given an iV x (complex) matrix A of rank A — 1, we 
would like to produce a nonzero vector in the nullspace of A, depending analytically 
on the entries of A. In general there is a topological obstruction to this; see Appendix 
\K\ However, the following result will be enough for our purposes. 

Fix N ^ 1. Let Mat(A^) be the space of all complex N x N matrices. We 
denote an A^ x A^ matrix by ^ e Mat(A^). We say that a map F : Mat(A^) 
is a polynomial map if the components of T{A) are polynomials in the entries of A. 
Polynomial maps are therefore analytic in the entries of A. 

In this section we prove 

Lemma 8.5. There exist polynomial maps Tjk ■ Mat(A^) — > C^, where j,k = 
1, . . . , A^, with the following property: 

Let A e Mat(A^) have rank N —I. Then all the vectors Tjk{A) belong to the nullspace 
of A, and at least one of these vectors is non-zero. 

Proof of Lemma \8. 51 ' We begin by setting up some notation. Given A e Mat(A), 
we write A^^''^^ to denote the matrix obtained from A by deleting row j and column 
k. We write Col(^, k) to denote the fc*'' column oi A. li v = {vi, . . . , vn)^ g is 
a column vector, then we write vj to denote the j*'' coordinate of v, and write {i^*'^ 
to denote the column vector obtained from v by deleting the A:*'' coordinate. Thus, 

From linear algebra, we recall that 

For any A e Mat{N) of rank A^ - 1, we have det (^A'^^ '^^^ ¥= for some (j, k) . 

(8.12) 

For any A e Mat(A^), and any v e Nullspace(j4), we have 

^(j,fe)£)(fe) = -[col{A,k)]^ ^^'> vk. (8.13) 



(Equation (|8.13p expresses the fact that {Av)i = 0, for all i ^ j.) 
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If rank(^) = iV - 1 and detA'^^'''^ 0, then the space of solutions of (ISTT^ . and 
the nuUspace of A are one-dimensional; hence, in this case the nullspace of A consists 
precisely of the solutions of (|8.13l) . 

We now define V j^^i^A) e to be the element, u, in the nullspace of A^ whose N 
components are constructed as follows: 

Set z;fc = rdetyl(^''=)j^ (8.14) 
and obtain the other TV — 1 entries comprising the vector v^^^ by solving 
^(j.fe) = - [col(A, fc)] ^ [detv4(J''=)]^ . (8.15) 

If det 7^ 0, we can solve (|8.15p uniquely for v^^^ by Cramer's rule and together 

with (|8.14p construct v. 

Note that each component of the vector Fj-fe (A) has the form 

det ^^'^^ X Polynomial in the entries of A . 

In particular, A i— > V jk{^A) is a polynomial map. Furthermore, Fjfc(^) = if 
det = 0. Also, V^k{A) if Ae,iA^^-^^ ^ 0, since the k^^ coordinate of T.jk{A) 

is [det^(J''=)]^ 

If rank(A) = TV — 1, then some Vjk{A) is non-zero, thanks to (|8.12p . Moreover, 
each Tjk{A) always belongs to the nullspace of A. Indeed, fix j, fc; if det^'-^^'^^ = 
then Tjk{A) = e Nullspace(A). If instead det A^^'''^ ^ 0, then Nullspace(A) consists 
of the solutions of (|8.13p : and we defined Fj^ to solve (|8.13p . Thus, in all cases we 
have Tjk{A) e Nullspace(A) if rank(v4) = A^— 1. The proof of Lemma [5751 is complete. 

8.4. Hamiltonians Depending on Parameters. Recall the operator H^'^^ = 
— A + £V^(x). In this section we complete the continuation argument (and the proof of 
Theorem lS.ip by showing how to continue Properties I. - IV. , listed at the beginning 
of section [8.11 beyond any critical value, £c, where one of these properties may fail. 
This argument is based on appropriate application of Lemma 18.41 and Lemma 18.51 



Let £c and /i^ be as in Definition! 
Without loss of generality we can assume K* = K. We work in the Hilbert spaces 

^kr = j E P"^"" + ^e^'"'^'" " + re^'K^ H : ^ 1^(^)1' < ool (8.16) 

2(l + |mn«|c(m)|2<cx)l (8.17) 



me 5 



We will apply the results of section \8.2\ and section \873[ with the analysis centered 
at (fic,£c) rather than at (0,0); see Remark \8.2l We shall use that 

H^^^ : ^ ^ is self-adjoint for e real . 
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Let M be a positive integer, chosen below to be sufficiently large. We regard ^ as 
the direct sum L^^ @ L^i, where Lf^ consists of Fourier series as in (|8.16p such that 
the c(m) = whenever |m| > M, and L^j consists of Fourier series as in (|8.16l) such 
that the c(m) = whenever |m| ^ M. Similarly, we regard ^ as the direct sum 
^lo ffi ^hi using (iSTTl) . We set N = dim{LfJ. 



Let Hio and Hhi be the projections that map a Fourier series as in (|8.16l) or (|8.17p 
to the truncated Fourier series obtained by setting all the c(m) with |m| > M, or 
with |m| ^ Af, respectively, equal to zero. We may view H^^^ as the mapping 

i^lo ) ^ [ C(-) D(^) ) [ ^lo 

with A(^) = Hhi H*-"^ Hhi, B*-^^ = Hhi H'-^^ Hio, 
C(^) = Hio H'^'^ Hhi, and D^"'' = Hio H^"^ Hio . 

By choosing the frequency cutoff, M, to be sufficiently large, we have 



/jjj Ljjj has a bounded inverse; say 

-1 



C . 

Therefore, for all in some fixed small neighborhood of {fic^c) we have 



^ C . 



1p\o 



The eigenvalue problem 

H^^'^xp = ^ ip for ip 
is equivalent to the system 

That is, 

^hi = - (a^'^ - fil) ' b(^) Vio, and 



hi 



(8.18) 



(8.19) 



.20) 



^hi, C(^)Vhi + i^^'V^lo = M V-lo . 



0, 



.21) 
.22) 



where we regard A^^") — as an operator from iJj^j to L^j. Note also that B^'^^ipio e 
since V'lo e if^ hence {A^^"^ - ^lY^ B^"^ V-'io e i/^j, thanks to (|8.19p . which holds 
under our assumption that (/i,e) is near (/ic,£c)- It follows that 



by the definition of C^^^ . 
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The operator in square brackets in (|8.22p will be denoted as V{^,£) . (8.23) 

Thus, 'D{ijL,e) is analytic in (/i,£), where (/i,e) varies over a small disc about (/ic,£c) 
in C'^. We may regard I?(/i,e) as an iV x matrix. Thus, 

/ _ / V'lo \ eigenfunction of i/^'^-* with eigenvalue a 

V ■0hi / 

if and only if (I8.2ip holds and i/^io is a non-trivial solution of V{iJ.,e)ip\o = 0. 

(8.24) 

It follows that fi is simple, i.e. a multiplicity one eigenvalue of H^'^^ if and only if the 
N X N matrix P(/i, e) has rank iV — 1. 

We shall now apply Lemma lSTSl to T>{fj,, e) e Mat(iV), the space oi N x N complex 
matrices. Let Tjk denote the polynomial map given by Lemma f8.5l For j,k = 1, . . . , N 
we define 

^|^{^i,e) = r,fe(2)(/i,e)) 

and set 

= '[a'^'^ B'^^^ijl^ 

as in dl^. 

By Lemma [8.51 l|8.2ip and (|8.22p we now know the following for 

/ W<f(M,£) \ 

V'^"'=(/i,e) = 

(A) il}^^{^,£) e -ffK.r depends analytically on (^,£), 

for (^,e) in a small neighborhood of {fic^c) (8.25) 



(J5) If is a simple eigenvalue of H^'^\ {fj. — /Zc, £ — £c small), 
then all ^^''{^, e) are in the nullspace of i/*-'^-* — /i/ . 
Furthermore, at least one of the ^■''^{fj,,£) is non-zero 
and is therefore an eig enfmiction of (8.26) 

Let us write out the Fourier expansions of the ip^'^{fi, e). We have 
[^p'''i^l,£)] (x) = c^'''^(m,/x,£) [ e'^"-" + fe'^^'"'' + re^^'^-.x j (g 27) 

me 5 

The coefhcients c^''{m, ii,£) depend analytically on {fj,,£) G U, where t7 is a small 
neighborhood of {fic,£c), which is independent of m. Moreover, since ^'■"^(/J,, £) is an 
analytic ^- valued function, it follows that 

(l + |m|^)^ |c-''''(m, ^, £)|^ is bounded as (^, £) varies over J7 . (8.28) 

meS 
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(Perhaps we must shrink U to achieve (|8.28l) .) 

With a view toward continuation of the Properties I. -IV. (enumerated at the start 
of sectionlS]) as e traverses any critical value (Definition [521), we state the foUowing 

Lemma 8.6. Suppose there exists a sequence of eigenvalues (/Lti/,ey) —>■ {hcEc) 
with < £y < Ec, such that for each v the following properties (Al )-(A4:) hold: 
(Al) is a simple eigenvalue of H^^'^^ on L\^^, with eigenfunction 

me 5 

(A2) III, is a simple eigenvalue of H^^"^ on ^kt' ^^t^ eigenfunction 

me 5 

(A3) £i{fj,iy,e^) ¥= 0, i.e. ii^, is not a eigenvalue of H^'^'^K 

(A4:) The following non-degeneracy condition (X^ ¥= 0) holds: 

2 wim) [c,(m)]' ^ 0, (8.29) 

me 5 

where {w(m)}me5 ci^s fixed weights, such that 

|w(m)| s£ C(l + |m|), me 5. (8.30) 
Our choice of weights (see (|4.ip ) is: 

Then, there exist a (non-empty) open interval I = (ec,Sc + S), a real-valued real- 
analytic function 13(e) defined on I , a function ip^ e ^ depending on the parameter 
eel, and a countable subset C c I, such that the following hold: 

(i) (-A + eVh) ip^"^ = /3(e)v?(^) for each e e Z. 

(a) hm_^^_^+ I3{e) = /^c- 

(Hi) C has no accumulation points in I, although Sc may be an accumulation point 
ofC. 

(iv) For each e in T\C 

(a) j3{e) is a simple eigenvalue of —A + eVu on ^, 

(b) 13(e) is not an eigenvalue of —A + eVh on ^, and 

(c) the quantity A| , arising from the eigenfunction tp^^'^ via formula (|4.ip 
(with <i>i replaced by (p*-'^-*) is non-zero. 

Proof of Lemma 1 8. 6\ RecaU that the zeros, /i, of the renormahzed determinant, 
£'i(^,e), defined in section FTSI are precisely the set of i eigenvalues of H^'^\ 
Thus, tracking the set of (/i,£) such that Assumptions {Al)-{A4:) and in particular 
(^3) suggests that we introduce, for (/i, e) e U, the matrix- valued function: 

F^''{li,e) = ( 2 wirn) [c^'' {m, fi,e)f ] x fi(^,£). (8.31) 

V me5 / 
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F^'^{lj,,e) is an analytic function on U. We define 

)j,fc = l,...,Af 



F(/i,e) = {F^'{^l,e) )^.,^^ ^ (8.32) 



Thus, F : U ^ is an analytic map. 

Now for each i/, (|8.26p applies to {fXy^e^) , since is a simple eigenvalue. Thus, 
for some j/c, the function 'ip^^{^jy,ey) is a non-zero null-vector of H^^"'^ — /x^/, i.e. 
an eigenfunction of H^^'^\ Since by hypothesis ^i, is an eigenfunction of H^'^"'^ sat- 
isfying (|8.29p with eigenvalue /i^ and since /i^ is a simple eigenvalue of H'^^^\ the 
corresponding eigenfunction Vl^y satisfies: 

'^u = lu '(p''''{fJ-u,£iy) for a complex constant 71, 7^ . 

Therefore, 

7^ 2 u'(m) [c^(m)]^ = 7^ 2 [c^'^im; n^,e^)f 

meS me 5 

= 7' ^'"/^'^'^^^ see dOJ) and (|H2S1). 

The second equality holds since £i{fi^,e^) 7^ 0; see hypothesis (^3). It follows that 
F^'^{liv,£u) for some jfc, i.e. 

F{fi^,e^) 0, for each 1/ , (8.33) 

where {e,^} is a sequence tending to £c from below. 



We complete the proof of Lemma 18.61 by application of Lemmata 18.41 and 18.51 
for appropriate choices of P(fi,e) and F(fj,,e). Let S-rifJ-jS), denote the renormalized 
determinant ([7^ . Let P{n,e) = fr(M,e) and F{fj.,e) be given by (jOTI) . ((O^ . We 
now check the hypotheses of Lemma 18.41 First note 



P(/i, e) = if and only if /i is an t eigenvalue of iJ^^-*, and the multiplicity of 
/X as a zero of P{fJ-, s) is equal to its multiplicity as an eig envalue of i/^^). (8.34) 

Because i?^^^ is self- adjoint for real e, we see from (|8.34p that 



if {fj.,e) e [/, £ is real, and P{fj.,e) = 0, then fieR. (8.35) 
Moreover, from (|8^ . ((833| and (|834)) we see that 

ifJ-v, £u) e U for each 5= 1, (/^i/, e^) ^ (/Xc, Ec), as ^ co; (8.36) 

and for each v, we have 

P{n^,e^) = 0, d^,P{fi^,e^) ^ 0, F{^^,e^) ^ 0, e R, < eq • (8.37) 

Recall that F : U ^ is an analytic mapping and P{n, e) : U ^ C is analytic. 

Results ((05|) - ((07| tell us that conditions (Al) and (A2) of section hold for 

our present choice of F and P. Therefore, Lemma 18.41 applies; see also the remark 
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immediately after its proof. Thus, we obtain a positive number S and a real-analytic 
function f3{e) such that the following hold: 

For each e e (ec,ec + <5), we have P(^(e),e) = and d^P{(3{e),e) 7^ 0. (8.38) 
Moreover, for all but countably many e e {Ec, ec + 5), 

with their only possible accumulation point at Sc, we have F(/3(e),e) 7^ 0. (8.39) 

By and ([QI)) . we have 

For each e in (£c,£c + <5), the number /3(£) is a simple r eigenvalue of i/^^-* . 

(8.40) 

Therefore, from ()8.26p we have that 

For each e in (cc, Ec + (5), all the ip^''^^, e) are in the nuUspace of i/^*^-* — l3{s)I. 

(8.41) 

Recalling (|8:27l) and ([OT|) . (|832l) . (lOOl) . we now see that 

For each e in (£c,£c + <5) outside a countable set (8.42) 
with its only possible accumulation point at Sc, £i(/3(e),e) ¥= 

and there exists jfc such that ^ u'(m) [c-''^(m, /3(e), e)]^ 7^ 0. (8.43) 



Now unfortunately the pair (j, k) in (|8.43p may depend on e. However, (|8.43p 
implies that for some fixed (j, k) = (jo, fco), the function 

^ w(m) [c'o'=n(m,/3(e),£)]^ (8.44) 

me 5 

defined for e e {Sc,£c + <5) is not identically zero. Since this function is analytic in e, 
it is equal to zero at most at countably many e. Moreover, the zeros of the function 
(|8.44[) in {ec, Sc + S) can accumulate only at Ec and at Ec + S. By taking S smaller, we 
may assume that the zeros of the function (|8.44p can only accumulate at Ec- 
We now set, for e e {ec, Ec + S): 

^(e) ^ ^'JO'feo(/3(e),e) e and 



<lr('^)(x) = cJ"'''"(m;/3(e),e) 



me 5 



We now have that the function /i = /3(e) satisfies Properties I. -IV. for all e G 
{Ec,Ec + 5) except possibly along a sequence of "bad" e's which tends to Ec- Properties 
I.-IIL, that /3(e) is a eigenvalue in each of the subspaces ^ and ^, and not an 

^— eigenvalue, hold for all e e {ec,Ec + S), except possibly along the above sequence 
of bad e's. This completes the proof of Lemma \8M 



To complete the proof of Theorem 15.11 for e of arbitrary size we require 

Lemma 8.7. For all e e (0, co) outside a countable closed set there exists a 
Floquet-Bloch eigenpair /i e M, (p e ^ for —A + eVu, with the following properties: 

a. ^ CqE + Ci, where Cq and Ci depend only on Vh- 

b. fi is a multiplicity one eigenvalue of —A + eVh on ^. 

c. /i is not an eigenvalue of —A + eVh on ^. 
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d. The quantity A|, arising from ip by formula (|4.ip is non-zero. 

Theorem 15. II is an immediate consequence of Lemma 18.71 which precludes option (2) 
of Lemma [Ol and Proposition 14. II 

Proof of Lemma \8. 7\ Set Co = max \Vh\. By our the analysis of section IH there 
exists e" > 0, a sufficiently large constant Ci, such that for e e (0, there exist /i, ip 
satisfying (a.)-(d.) . 

Now suppose that Lemma [5771 fails. Then, by Lemma I5TT] there exists Sc e (0,oo) 
such that for all e e (0,ec) outside a countable closed set, there exist fi,(p satisfying 
(a.)-(d.) but 

for all e;", > Ec, assertions (a.-d.) fail on a subset of {0,s].) 

that is not contained in any countable closed set. (8.45) 

We will deduce a contradiction, from which we conclude Lemma 18.71 

By assumption, we can find a sequence ei < 62 < ■ ■ • < < ■ ■ converging to Sc, 
such that each e„ gives rise to a Floquet-Bloch eigenpair fi^, e K, ip^ e ^ satisfying 
properties (a.-d.) . Thanks to (a), we may pass to a subsequence, and assume that 
A*!/ ^ A*c as f ^ CO, for some real number /ic, with 

|a*c| Co Ec + Ci (8.46) 

Since the Floquet-Bloch pairs /Xi/, ipi, satisfy (a.-d.), and since e^, \ Sc and ^ ^ 
as — > 00, Lemma 18.61 applies. Thus we obtain a non-empty open interval I = 
{Sc,Sc + S), a real-valued real-analytic function /3(e) defined on I, a function ip^^'' e 
-^K T parametrized by e e I, and a countable closed subset C c I satisfying properties 
(i)-(iv) of Lemma [8.61 We will prove that 

\l3{e)\ s; Co £ + Ci, for ah e e Z = (e^ £c + (5). (8.47) 

Once (|8.47p is established, we conclude that we can satisfy all assertions (a.-d.) of 
Lemma [5771 for all e e T\C, by taking /i = /?(e) and ip = (p^'^\ However this contradicts 
property (|8.45p of Sc- Thus, it suffices to prove the bound (|8.47p . 

To establish ((5T7)) . we fix e e I\C = (ec,£c + ^)\C. For any e > 0, let Ai(e) ^ 
A2(e) ^ . . . denote the eigenvalues (multiplicity counted) of —A + eVh on ^. Then, 
since /3(e) is a simple eigenvalue ( assertion (iv) of Lemma 18.61 ). there exists k such 
that 

^(e) = ^fe(^) < Afc+i(e), and Afe_i(£) < Afe(e) unless = 1. 
Fix ^ > such that 

hi^) < h+ii^) - V, and Afc_i(e) < Afe(e) - fj unless k = 1. 



From the min-max characterization of eigenvalues, we have the Lipschitz bound 
|Afe(e)-Afe(e')| |e - e'| • max |14| = Co-|£-£'|, (8.48) 
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for any e, e' > and any k ^ I. Hence, as e varies in a small neighborhood of £, we 
have 

\X-,ie) - X-,ie)\ ^ Co-\e-e\, (8.49) 

and also 

Afc+i(e) > Afc(£) + if/, and X^^iie) < Afc(e) - ^fj unless fc = 1. (8.50) 

We have taken e e T\C. As e varies in a small neighborhood of e, we have e e 2r\C, 
thanks to property (iii) of Lemma [8.61 Therefore, /3(e) is an eigenvalue of —A + eVh 
on ^, i.e. I3{e) = Afc(e) for some k. We now show that this value must be fc. 
Since /3(e) is a real-analytic function of e, and since /3(e) = Aj.(e), we have 

Afc(e-)-i?7 < /3(e) < Afc(e-) + ^77 (8.51) 
for all e close enough to e. From (I8.50p and (|8.5ip , we have 

Afc_i(£) < ^(e) < Afc+i(e) 

and therefore, /3(e) = Aj.(e) for all e close enough to e. Estimate (18.49^ now shows 
that the real- analytic function /3(e) satisfies 

< Co, for e = e. 

Since e was taken to be an arbitrary point of T\C, and since T\C is dense in I by 
(iii) of Lemma 18.61 we have 

< Co, for aU eel. (8.52) 

Recall that I = {Ec, Sc + S). Our desired estimate (|8.47p now follows at once from 
(|8.46p . (|8.52p and (ii) of Lemma l8.6l The proof of Lemma l8.7l and therefore of Theorem 
15. H is now complete. 

9. Deformed honeycomb structures and robust conical singularities. 

In the previous sections we established the existence of conical singularities, Dirac 
points, in the dispersion surface for honeycomb lattice potentials. These Dirac points 
are at the vertices of the Brillouin zone, Bh. In this section we explore the structural 
stability question of whether such Dirac points persist under small, even and A^- 
periodic perturbations of a base honeycomb lattice potential. We prove the following 

Theorem 9.1. Let V(x) denote a honeycomb lattice potential in the sense of 
Definition \2.1\ Let M^(x) denote a smooth, even and Ah- periodic function, which 
does not necessarily have honeycomb structure .symmetry, i.e. Wijx.) is not necessarily 
TZ- invariant. Consider the operator 

H{r]) = -A + F(x) + r]W{x) , (9.1) 



dm 

de 



dm 

de 



* For A: = 1, we have /3(e) < A2(e). 
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where rj is a real parameter. Let k = K, be a vertex oj Bh- Assume that for rj = 0, the 
operator H{0) has an - eigenvalue, /i(K*), of multiplicity two, with corresponding 

orthonormal basis {$i,$2} with $i 6 ^ and $2(x) = <i>i(— x). Assume Aj, given 
in (|4.ip . is non-zero. Then, the following hold: 

1. There exist a positive number rji and a smooth mappings 

,y ^ ^i^) = ;i(KO + 0{ri) eR andf]^ K^''^ = + 0(ry) e Bh, (9.2) 
77K^05"^(x;K('')) = <f>,{^) + Oiv)eL^iR/Ah) (9.3) 

defined for \r]\ < rji, such that H(t]) has an i^(^) eigenvalue, /i'-''-*; of geomet- 
ric multiplicity two, with corresponding eigenspace spanned by 

2. The operator H(ri) has conical-type dispersion surfaces in a neighborhood of 
points kI''-* = K* + 0{rj) with associated band dispersion functions, /Lt+''(k), 
defined for k near kI'''* ; 

^V'''(k) -^(kI")) = 77b(") ■ (k-Kl"') + (^Q''')(k-Ki")))' 1 + (k-Ki"') 

(9.4) 

^^'''(k) -M(Ki''') = r^b*"' ■ (k-Ki")) - (^Q('''(k-Ki"')) ' 1 + si"' (k - kI"') 

(9.5) 

where 

• b*^''^ 6 depends smoothly on rj. 

• Q^^\-) is a quadratic form in k = (ki, K2) e M^, depending smoothly on 
rj and such that 

{\h?-CU){nl + Kl) < Q{n-^) < (|A„p + qr;|)(.c? + «;2) (9.6) 

/or \ri\ < 771 and k = (ki, K2) e R^, with rji small, and 

• E^\k) , E^^\k) ^ C \k\ for \r]\ ^ 771 anc? |k| ^ where k (small) 
and C < CO are constants. 

Remark 19.21 below shows that Dirac points are unstable to typical perturbations, 
W e C=^(]R7A/i), which are not even. 



Remark 9.1. For ri = 0, Theorem \9. 1\ reduces to Theorem \5.1\ which covers the 
case of the undeformed honeycomb lattice potential. In particular, if the perturbation 
W is itself a honeycomb lattice potential, then kI''-* = K, . 

We now prove Theorem l9.1l As earlier, without loss of generality, we assume K* = K. 
The family of Floquet-Bloch eigenvalue problems, parametrized by k e Bh, is given 
by: 

iJ(k;?7) 0('')(x;k) = ^(''^(k) (x; k), where 
iJ(k;77) = -{V + ikf + V{^) + 7] VF(x) . (9.7) 
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By hypothesis, we have: 

For 7y = 0, H{K; 0) = — (V + iK.)^ + V{x.) has a degenerate eigenvalue n{K.), 
of multiphcity two with 2-diniensional L^- eigenspace: span{(/)i(x), (/)2(x)}. (9.8) 

Introduce the projection operators 



Q\\f = l]<'^„/>0j(x) andQj. = /-Oil . (9.9) 
i=i 

We seek solutions of the Floquet-Bloch eigenvalue problem in the form: 

2 

,/)('')(x;K") = 2«.<^.-W + V(t>^''H^), Oii^^^'") = , (9.10) 

K^") = K + ryKi'" , (9.11) 

^(v) = ^(0) + ^^ii..v) ^ (9.12) 

dim null space H{rj) - fi'-'''h ) = 2 . (9.13) 
Substituting these expansions into (|9.7p yields: 

i7(K)-^(°)/) 0(1^") + ( -2zKi^" • Vk + - ^i^^'"^ - 7?|Ki''f ) 0(1'") 

= - ( -2zKi'" • Vk + - Ai^i'") - JylK^'"!') 1^ 2 a,^, ^ . (9.14) 
We take the inner product of (|9.14|) with 0i and 02- This yields the system: 

+ ?7 (<^i,(-2jK(^'"' ■ VK + W-jQi.^'''"') = , (9.15) 

+ »7 (02, ('-2jK<i'''' ■ Vk + W-) Qi.^'^'"') = . (9.16) 



In obtaining the system (|9.15p - (|9.16p we have used 
1. 

Vk = (Vx + iK) = e-'^'Vxe''^", (9.17) 

2. W is real-valued and thus —2iK^'^ • V + is self-adjoint, 

3. -2i<0i, K^-" • Vk 02> = -2i<$i, K^-" ■ Vx -I>2> = (k'"''' + iK*"'''' 
where = e'^"0j, by Propositior 

4. 0110(1''') = and 0x0(1''') ^ ^(i,r,) 



An equation for 0(1''') is derived by applying Ox to (|9.14p and using 
Ox0j = 0, j = 1,2 and 0x0^'^"^ = 0^''"^- 
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= -Qi_ ( -2iK'^'' • Vk + H/) yj] j ■ (9-18) 
Introduce the projections 

2 

gi|F= 2<'i'j,^><fj, Q± = i-Q\\ (9.19) 

and note that these projections satisfy the commutation relation 

e*K x Q = Q e^K-x 20) 

Using (|9.17p we rewrite (|9.18p as an equation for 

2 

/:(M,Ki^",ry) = - ^ Qx ( -2iKi^" • Vx + H/) where (9.22) 

£(m,K1'",77) = -A+^-^(") +77 Qx (-2^Kl■"• Vx + Ty'lKi'^l^) Qx (9.23) 

Note that C is self-adjoint on Q±L^. For 77 sufhciently small, we have that C : 
Q±H^ Q±Ll^ is invertible. Solving yields 

2 

$(1,^) = a, C{fi,K''\r^y'Qi_ { ~2iK''^ ■ + W) $j . (9.24) 

Using ((9?2T|) to express ^^^'''^ in terms of and substituting (|9^ into (|9?T5|) - 

()9.16p . we obtain the following linear homogeneous system for aj, j = 1, 2: 

^ (""■•"■«■■"■")(:) - (0) ■ •''■^^) 

where 

/ - (^(i.'J) +,7|Ki'''|2) +<$i,Ty$i> + r? an (iTj''' + 1X3''') + <*i, Ty*2> + \ 



y Aj (a-^^'" + jXj^'") +<<t.i,VK'I'2> + »; 6 - (/xti'") +ij|Ki'''l2) +<$i,W$i> + ,, a22 y 

(9.26) 

ail, 122 and b are functions of fi^^'^\ K.^'^ and 77 and are given by the expressions: 

au = (g± ( -2iKi'" • Vx + ) $i , Qx ( -27K1'" • Vx + ) , ^ = 1, 2 

(9.27) 

b= (gx ( -2jk1'" • Vx + ) $1 , gx ( -2iK''^ . Vx + ) $2) (9.28) 
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where L = >C(/i, K^,?]) is defined in (|9.23p . Note that an and 022 are real. The 
matrix M. (/i*-^-*, K^, r/) has the structure: 



(9.29) 



where A;; and B are smooth functions of [jS^^ ^Yi}'"^ ^r\) , which can be read off from 
(19:261) and (lOTll - dOSl) : 

= <$i,VK$,>-?7|Ki|2 + r;a„, ? = 1,2 (9.30) 
B = T^{Kl + iK^) + <$i,VK$2> + ?? 6 . (9.31) 

A consequence of the above discussion is 
Proposition 9.2. 

1. The pair {11^"^^ (j)^^"^) is anL^(,,)- eigenpair of H{r]) = —A + V + rjW, where 
f] is real and K^'') e R^, ^(''), e i^,,, are defined in (|9J0l) - (|932l) . if and 
only if 

detM (^i'^^\K\ri] = . (9.32) 



2. By self-adjointness, for ry e R and 6 M^, z/ /x^^-* is a solution of (j9.32p 

3. pS^^ is a geometric multiplicity two - eigenvalue of H{jf) if and only if 
the triple (/i^"'^^ K"'^, ry) is such that the 2x2 Hermitian matrix, M ,ri^ 
has zero as a double eigenvalue, i.e. M. (/i^"'^-*, K^, 77) is the zero matrix. 

Now, up to this point we have not used the hypothesis that W^(x) is an even 
function (inversion symmetry). We now impose this condition on W. For the case 
where W is not even, see Remark 19.21 at this end of this section. 



Claim 1: l^(x) = VF(-x) =^ <$i,T^$i> = <$2,H/$2> and an = 022. By 
(lOO)) . it follows that An = A22. 



Proof of Claim 1: Recall, as in Theorem 15. 11 that $2(x;K) = $i(— x;K) and thus 
{^2,W^2) = j $i(-x)iy(x)$i(-x)dx = j |$i(x)|2T^(-x)dx = <$l,W^$l>. 



Furthermore, one checks easily that an = a22. 
In this case, we set 



an = a22 = a = a(/i(i),Ki,77) and An = ^22 = ^ = A K^, 77) . (9.33) 

Here a = an = 022 and b are functions of p^^\K.^ ,ri, displayed in (|9.27l) and (I9.28p . 
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By Proposition 19.21 and the above Claim 1, if T^(x) = W{—x.), then we obtain a 
double eigenvalue if and only if 

H^^^ ~ A(^n'^^\K\f]'^ =0 and S (^^^^ K\ 77) = 0. (9.34) 

By analyzing the solution set of (|9.34l) for small 77, we shall prove the following: 

Proposition 9.3. For each realrj in some small neighborhood of zero, there exists 
a unique = K^''' = (^K^'' , K^'"^^ and /^(i) = /i^^^'') such that /i^'') = ^(0) + jy/i^i''') 
(see (|9.12p ) is a geometric multiplicity two L^(R^/Ah)- eigenvalue of H{rj]\i}'^). 
Proof of Provosition \9.Sl - Consider (|9.34p for fi^^'"^^ and K^'** for 77 = 0. We have 

/i(i'°)-^(Ai^i'°),K\o) =0 ^ = <$i,W^$i> and (9.35) 

B (/x(i^°),K\o) = ^ Ti (^kI° + iK^/'^ = -<$i,W^$2>- (9.36) 
Equation (|9.36p is equivalent to the two equations: 

i^i." = -SR ( ( a; )-i <$!, W^2) ) , if2'° = -3 ( ( a; )-i <$i, W^$2> ) (9.37) 
We next consider the case 7; 7^ 0, real and sufficiently small. 

Claim 2: A and B, defined via dSSSJ-dSSHl), are smooth functions of (/z^^^ , 77, K}'°). 
Moreover, there exist constants ci > 0, di > and 770 > such that for all |K^ — 
K^'"| < ci, l^(^) - ^(^'"^l < di and I77I < 770, we have 

1. A = (^,,W^,} + 77/^(/i(i),Ki,77) 

2. B = Yi (Kl + iK^) + <$i,M^$2> + V fB{ti^^\K\7j) 

3. /a, fB = 0(1) . 

4. = 0(r,), \8f,A\ < 1/2. 
We leave the verification to the reader. 

An immediate consequence of Claims 1 and 2 is: 

Claim 3: Assume Aj ^ 0. Then, for \r)\ < 770, \K^-K^'°\ < a and - ^(1^°) | < di 
equations (I9.34p are equivalent to the system: 

Ki'^ + = (A;)"\$i,W-<i>i> + 77^^(1), 77, 7^«,i^«), (9.38) 

= <$i,l^$i> + 77/A (/i(i),Ki,77) . (9.39) 

By the above, we have 

Claim 3: //("J) = + t^^^^''') is an L'^ {M.'^ / A,,)- eigenvalue of iJ(K(''\ 77) of geometric 
muhiphcity two (see (|97[0l) - ([9l^ ) if and only if ^^i-''), if^''' and KI'"^ satisfy: 

X^-" - 5R( (AJ)"\$l,W^<i>l> + 77^^(A^(l■''^77,i^^^i^2'''')) = 
^2'" - »( (AJ)"\$l,l¥<i>l> + r;i^(A^(^■''),77,X^^if2'''')) = 
^(M) _ (<$i,W^<i>i> + r, fA{pi^'''^\v,Kl-\K'^^^)) = (9.40) 
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So in order to prove Theorem 19. II we seek a solution of (|9.40l) in a neighborhood 
of its solution for r/ = 0, given by (j9.35p and (|9.37l) : 

Note that the right hand side of ()9.40p defines a smooth map from a neighborhood 
of {h^^'°\kI''^,kI'^) to R3 with Jacobian at dMB, for = 0, equal to the identity. 
Hence, by the implicit function theorem, there exist a positive number, 771, and smooth 
functions: 

ry /i^i'"), T^^K^''^ = (kI'\kI''''^ , (9.42) 

defined for I77I < 771, such that /i^^'''^ K^'** is the unique solution of (I9.40p for all 
1 77 1 < r]i in an open set about the point (|9.4ip . This completes the proof of part (1) 
of Theorem O 

To prove part (2) of Theorem 19. 11 we need to display a conical singularity in the 
dispersion surface about the point (K'^''^ , /i'^'') ) , (|9.2p . For this we make strong use 
of the calculations in the proof of Theorem 14.11 In particular, —A + 1/, /i^"-*, K and 
4>j, j = 1,2 of the proof of Theorem 14.11 are replaced by H{r]), ii^^^K.^"^^ from the 

proof of part (1) and {0^''"' , ^j''"'}, now denotes an orthonormal spanning set for the 
i2(R2/Aft) nullspaceofir(7?;K(''))-^('')/. Then, {^["'\ <i>'^^^ } = {e*K-^0i''\ e'^'^^''^} 
is an orthonormal spanning set for the i^(^) nuUspace of H{r]) — /i^''^/. Note also that 

= the Floquet-Bloch states associated with the unperturbed honeycomb 
lattice potential, V. 

We must study the Floquet-Bloch eigenvalue problem (compare with ()4.4p - (|4.5p ) 

- (Vx + i (k('') + V{x) + r7VK(x)j V'^'')(x;K('') + «;) 

= /i(K(") -I- k) V^"^(x; K^") + k) , (9.43) 
v;K('') + k) = V('')(x;K('') + k), for all v 6 A . (9.44) 

We express = /^(''^(K^'') + k) and V^''^^; K^*?) 4- k) as: 

^ = ^('') + ^(1''')(k), 
2 

^W(x;K('') + At) = 2 aj0?'\x) + V^^i'^Hx), (9.45) 

where (K^**)) denotes the perturbed double-eigenvalue constructed above 

with corresponding orthonormal eigenfunctions 4>[^\'x), j = 1,2. 

Precisely along the lines of the derivation of (|4.19l) in the proof of Theorem I4.1i 
we now find that for \k\ small ^('')(K('') + k) is an eigenvalue of the spectral problem 
if Ai^^'''^ = A«(^'''H'«) solves 

detX(/x(i'''\K;?7) = 0, (9.46) 

where 

M{fi'^^^^\K;f]) = Mo{fi^^'''\i^;v) + Xi(m^''''\«;^), (9.47) 
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where 
and 



where M^j,{ri) are smooth complex- valued functions of r/. Note that 



(9.48) 



?7=0 



-Aj (ki + ZK2) 



^n(^), are real and Af2"i(77) = MUv) 

and furthermore 

Thanks to (|9.48l) . the equation dot A1o('^i ^^w) =0 is equivalent to 

+ { [mUv) + MUij)] • « } + dct [{M%{ii) ■ «)^. , 
The solutions have the form 



(9.49) 



(9.50) 



(9.51) 



where Q('')(k) is a quadratic form in k with coefficients depending smoothly on rj. 
For ?7 = 0, (j9.49p shows that the quadratic equation (j9.50p takes the form: 

u^-\\\^{Ki + 4)=Q; 

hence in (|9.5ip we have 

[ M°i(7y) + M^ri) ]\, = and Q^^\n) 

I ri=0 

Therefore, for |?7| < r/i (small) (|9.5ip takes the form 



+ (9.52) 



(9.53) 



where b'^''^ 6 depends smoothly on ?/, Q^^^\k) is a quadratic form in k (depending 
smoothly on ij) and 

W-C\v\){^i + 4) < 2^"^'*) < {\M\' + C\v\){ni + 4), (9.54) 

for I?/! ^ 771 and k = (ki,K2) e K^. Thus, the solutions of det A4o{i', k;?]) = are 
given by ((57551) and 
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We may now pass from solutions of det A4q{i>, k; t]) = to solutions of 
det M-{i^, k; r/) = as in the proof of Proposition 14.31 . Thus the -^k('!)+k" eigenvalues 
of H'^^^ are given by 



where E+{K;rj) ^ C\n\ for \r]\ ^ t^i, |k| k. The proof of Theorem 19. II is complete. 

Remark 9.2 (Instability of the Dirac Point and smooth dispersion surfaces). 
We here note a class of perturbing potentials, W , such that although —A + V has 
Dirac (conical) points, the operator —A + V + rfW has a locally smooth dispersion 
surface near the vertices of Bh- Assume that V is a honeycomb lattice potential, 
which is inversion- symmetric with respect to x = 0, i.e. Xq = m Definition \2.1l 
i.e. y(— x) = V(x). Let W e C*'(M), Ah- periodic, but without the requirement that 
W{x.) = W{-x) for all x. Then, typically (^i,W^i) <$2, VK$2>. In this case, 

see (|9.30p . For fi^^^ = fi^*^^ + rjfi^'^ to be an i^(^) eigenvalue, we found that it is 
necessary and sufficient that: 

or equivalently 



Thus, our eigenvalue equation becomes: 



A11+A22 , /Mil -^22^^ 



(1,.) ^ ^11 ^ ^22 ^ (g gg) 



When, A11—A22 0, each sign in (j9.56l) gives rise to an equation to which we may ap- 
ply the implicit function theorem to obtain a smooth function (K^,r/) 1— > ^(^■''^(K^, 77). 
In particular, at rj = 0, equation (|9.56p gives 



2 



= i^j |$x(x)|^ W{x)~W{~x) + I A^(j^i.o+,;^i.o)+<$,,Vy$,>|V (9.57) 

Therefore, for small rj, the two signs in (|9.56p give rise to two distinct solutions fi+'^^ . 
Thus, for small non-zero rj, the double eigenvalue disappears and the dispersion surface 
is smooth. 



Appendix A. Topological obstruction. 



Charles L. FefFerman and Michael I. Weinstein 
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In section [S] there arises the situation of an iV x complex matrix A varying 
within the space of rank iV — 1 matrices. It was of interest to know whether one 
can construct a non-zero nuUvector which is an analytic function of the entries of A. 
In this section we provide a 2 x 2 matrix counterexample that exhibits a topological 
obstruction. 

Let M cz Mat(2) denote the space of 2 x 2 complex matrices of rank 1. We prove 
the following 

Proposition A.l. There is no continuous map : M ^ C^\{0} such that 
(l){A) e Nullspacef^J for each AeM. 

Proof. Let cj) denote such a map. We proceed to derive a contradiction. For 

vectors w = ^ ^ e C^\{0}, define the 2x2 complex rank 1 matrix: 

A{v) = v®Jv = v{Jvf , (A.l) 

where J is skew symmetric and non-singular. Note: v i— > A(v) is a continuous map 
from C^\{0} to M. By skew-symmetry of J, A{y)v = and therefore 

Nullspace(A(w)) = C x w for each v e C^\{0} . 

Hence, for each v 6 C^\{0} there is one and only one non-zero complex number X{v) 
such that 

(l){A{v)) = X{v)v . (A.2) 

Since (j> is assumed continuous, the map v i— > \(v) is continuous from C^\{0} to C\{0}. 
Moreover, for aU v e C^MO} and BeR: A(e*%) = A{v). Hence, by (IX2|) 

A(e*%)-e'% = (/)(^(e*%)) = (t){A{v)) = \{v)v 

and therefore A (e*%) = e"'" X{v) . (A.3) 



Now for every 9 e = R/2ttZ and t e [0, 1], let 

v{e;t) = t ei + (1 - t) 63, where ei=(^J^, e2=(^J^ . (A.4) 

Note v{9;t) e C\{0} and introduce, for {9,t) e x [0, 1] 

ae;t) = X{v{e;t)). (A.5) 

We think of 6* i-^- (^{6; t) as a 1-parameter {t e [0, 1]) family of closed curves in C^\{0}. 

Taking i = 0, we have C{&; 0) = A (62) , for aU e S'^; and 

taking t = 1, we have C(6'; 1) = A (e^^ei) = e^'^X (ei) , for aU 9 e , 

by (jA.3[) . Thus by varying t between and 1 we obtain a continuous deformation of 
the unit circle to a point, remaining in C\{0}. This is impossible. □ 
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